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Type-set by plain 

Properties of neutral and charged anyon fluids are examined, with the main focus on 
the question whether or not a charged anyon fluid exhibits a superconductivity at zero 
and finite temperature. Quantum mechanics of anyon fluids is precisely described by 
Chern-Simons gauge theory. The random phase approximation (RPA), the linearized self- 
consistent field method (SCF), and the hydrodynamic approach employed in the early 
analysis of anyon fluids are all equivalent. Relations and differences between neutral and 
charged anyon fluids are discussed. It is necessary to go beyond RPA and the linearized 
SCF, and possively beyond the Hartree-Fock approximation, to correctly describe various 
phenomena such as the flux quantization, vortex formation, and phase transition. 
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1. Introduction 

Is a charged anyon fluid a superconductor? Is there any difference in its be- 
haviour from traditional superconductors described by Ginzburg-Landau-BCS the- 
ory? Can newly discovered high Tc superconductors be anyon superconductors? 
These are the main questions addressed in this article. Since Laughlin suggested 
that a high Tc superconductor may be viewed as an anyon fluid, an extensive study 
has been conducted by many authors by various methods. 

After three years of investigation, we can now have a coherent assessment of the 
current understanding. A fair statement, at the moment, is that a charged anyon 
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fluid seems to behave like a superconductor, but its properties have not been under- 
stood very well. We are not even sure if we have found a good approximate ground 
state of a charged anyon fluid, which should serve as an alternative to the BCS 
ground state of ordinary superconductors. The approximate ground state employed 
in the early investigation, characterized by a state of completely flUed Landau levels, 
might miss many of important phenomena such as the flux quantization, Josephson 
effect, and vortex formation. 

Anyons^^ exist in Nature. Excited states (quasi particles) in fractional quan- 
tum Hall effect (FQHE) obey fractional statistics. Laughlin's theory of FQHE 
implies that for a filling factor v = p/q, where p and q are coprime numbers, quasi 
particles have a statistics phase ±7r/g. The theory predicts a hierarchy structure 
in the Hall conductivity a^y as the filling factor or magnetic field varies, which has 
been confirmed experimentally at multiple levels. '^^"^^ 

Where else are anyons? Three years ago Laughlin made a bold hypothesis that 
anyons are in newly discovered high Tc superconductors.^^ 

There are three issues involved here. (1) First one has to show that anyons 
really emerge as excitations, starting from some spin, or electron, systems. There 
are arguments that the ground state of the Hubbard or t-J model in two dimensions 
near the half-filling may be in the flux phase or more specifically in chiral spin 
state, in which excitations arc anyons. So far only consistency arguments have been 
provided. (2) Secondly, assuming that there exist anyon fluids with a flnite density, 
one needs to know if such fluids exhibit a superfluidity or superconductivity. One 
has to midcrstand physics of anyon fluids. (3) Thirdly, one has to deduce physical 
properties of anyon fluids which can be subject to experimental tests. Specifically, 
one has to know whether or not (a part of) high Tc superconductors are anyon 
superconductors. 

The first issue is most difficult, and is outside the topics covered in this article. 
The third one is controversial, in both experimental and theoretical viewpoints. 
This article exclusively deals with the second issue. 

First Fetter, Hanna, and Laughlin showed, ^^'^^ computating response functions 
in the random phase approximation (RPA), that a neutral anyon fluid has a massless 
excitation, and that with electromagnetic interactions coupled the system becomes 
superconducting. Several months later, various authors confirmed the superconduc- 
tivity at T = by many different methods. Chen, Wilczck. Wittcn, and Halperin 
generalized the RPA analysis, deriving many physical consequences.^^ Wen and Zee, 
taking the hydrodynamic approach, gave a physical meaning of a massless excita- 
tion as a breathing mode of density wave.^° Canright, Girvin, and Brass's numerical 
analysis yielded a desired flux dependence in superconductors.'^^ The self-consistent 
field analysis of a charged anyon fluid by Hosotani and Chakravarty led to a new 
equation replacing the London equation in the BCS theory. Fradkin gave an anal- 
ysis of the model on a lattice. ^'^ 

At first all these analyses looked quite different from each other, though reaching 
to the same conclusion, namely a superconductivity at T = 0. In the mean time 
the investigation has been extended in various directions, which includes analyses of 
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finite temperature, vortices, conductivity, P and T violation, higher order radiative 
corrections, and so force. ^"^"^^ 

It has been recently shown^** that the RPA, self-consistent field method (SCF), 
and hydrodynamic approach yield the same results for many physical quantities such 
as the excitation spectrum and response function. In this paper we shall strengthen 
the statement. Wc show that RPA is exactly the same as the linearlized version of 
SCF, and that the hydrodynamic approach describes the same physics in terms of 
the density and velocity field as SCF does in terms of the gauge fields. Hence all 
these three arc equivalent. 

Quantum mechanics of anyon systems is precisely described in terms of Chern- 
Simons gauge theory. The essence of anyon dynamics is contained in the Aharanov- 
Bohm effect with respect to Chern-Simons gauge fields. We shall establish the 
equivalence between the two descriptions in the following three sections. It will be 
seen that the language of Chern-Simons gauge theory facilitates and simplifies all 
the discussions of anyon fiuids. 



2. Anyons 

Under the interchange of two identical particles, the Schrodinger wave function 
in quantum mechanics acquires a factor of either +1 or —1, depending on whether 
the particles are bosons or fermions. In two spacial dimensions there can be other 
possibility.^^ The interchange of two particles, say a and b, defines two paths Ci 
and C2 along which the particles a and b are transported to the original locations 
of b and a, respectively. Ci and C2 together form an oriented closed loop. Pick the 
paths such that none of the other particles are inside the closed loop. Depending 
on whether the loop is oriented in a clockwise or counter-clockwise direction, the 
operation defines P_(o, 6) or P+(a, 6). (Fig. 1) 

In three dimensional space there can be no distinction between P+(a, 6) and 
P-{a,b) {= P{a,b)), since P+ can be continuously deformed to P_. Hence 
P{a,b)^ = P+{a,b)P-{a,b) = 1, and P{a,b) = ±1. In two spatial dimensions, 
however, one can have 

P±(a, b) $(1, ■■■,q) = -e±'^^ $(1, ■■■,q) (2.1) 

where $(1, ■ ■ ■ ,q) is the Schrodinger wave function for a g-particle system. The 
minus sign in (2.1), retained for later convenience, may be absorbed in the definition 
of the statistical phase 6s- 6s — or -k {mod 2tt) corresponds to fermions or bosons, 
respectively. Otherwise the statistics of the particles is in between. It is said that 
particles obey fractional statistics. Such particles are generically called anyons. It 
is easy to see that (2.1) satisfies, for instance, an operational identity (Fig. 2) 

P.{a,b)P+{a,c)P+{b,c) = P+{a,c) . (2.2) 

If the wave function is to be well-defined in the limit two of the coordinates Xa 
and Xf, coincide, the identity (2.1) leads to 

$(l,---,g) =0 for -e*^' 1 . (2.3) 

Xa=Xi, 
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In other words, unless particles are bosons, the wave function must vanish when 
two coordinates coincide. Pauli's exclusion principle applies to anyons. 



Fig. 1 Interchange of two identical particles. The closed contour formed by 
Ci and C2 should not encircle any other particles. 



Fig. 2 The identity (2.2). 



A system of "free" anyons is described by the equation 



dt 

a=l 

where Va = d/d'x.a- The equation (2.4) must be solved with the boundary condition 
(2.1). This system defines a neutral any on fluid. We shall see below that except for 
the cases of bosons and fcrmions a neutral anyon fluid is not "free" . The energy of 
a many-anyon system is not the sum of single particle energies. An interaction is 
hidden in the nontrivial boundary condition (2.1). 

It is most instructive to go over to a new gauge. We define 



;(xi,- • • ,x„) = — tan 

^ TT Xal- Xbi 

a<b 



-1 Xa2 - Xb2 (2.5) 
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In terms of the new wave function, the equation and boundary condition become 

*|*' = i:-|;(v.-^B")(W))%' 

a— 1 

(Xa-X,)2 (2.6) 

= --Ve^'=a,« ln|x„-x,| 

P±{a,b)^'{l,---,q) = -^'{l,---,q) 

In the new representation the particles behave as fermions, but with a specific long 
range interaction described by B. is the Schrodinger wave function in the fermion 

representation. 

The B*^°)({xh}) term in (2.6) gives rise to two- and three-body interactions. It 
involves a velocity dependent potential. The interaction, which account for the 
anyon nature of the particles, can be interpreted as an Aharonov-Bohm effect or 
as a Chern-Simons gauge interaction, as we shall show in the subsequent sections. 
The gauge transformation potential fusing which connects the original $ and new 

is singular at = Xf,. One advantage of working in the new gauge is that the 
wave function is a regular, single-valued function of the coordinates {xo}. 

3. Aharonov-Bohm EfTect 

Suppose that there is a solenoid parallel to the a;3-axis at a;i = X2 = with a 
total magnetic flux A* = / dx\dx2 B^. Outside the solenoid there results a vector 
potential 

A^ = ^ , ^^ = , A^ = Q (3.1) 
2-Kr 

in cylindrical coordinates (a;i=rcos^, X2=rsva.(j)) or 

As far as E = B = 0, the motion of electrons outside the solenoid is not affected by 
the presence of the flux /x ^ in classical theory. 

In quantum theory a non- vanishing gauge potential A^, which locally generates 
vanishing field strengths, but is not globally a pure gauge, affects the motion of 
electrons. Let us focus on the two-dimensional motion of electrons in the X1-X2 
plane, supposing that a momentum in the a;3-direction is zero. The Schrodinger 
equation outside the solenoid (r > R) is 



_^ri 

2m \ r 



d-Zo-r + Ad-^'^) ho = i^*o . (3.3) 
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The general solution is 

oo 

$0 = ^ e''^ {aiJi-a{kr) + biJa-i{kr)} (3.4) 

l=—oo 

where a = e/i/27r?ic, fc^ = 2mE/Ti^ , and is a Bcsscl function of fractional order 
v. The wave function is supposed to vanish at the boundary of the solenoid. The 
wave function and energy eigenvalue depend on a, or on the magnetic flux /x in a 
periodic fashion. It is called the Aharonov-Bohm effect.''^ 

A basic assumption in deriving (3.4) is that the wave function $0 is a singlevalued 
function of r and (j). Let us define a new wave function by ^^{r, <p) — e*"'^<i?i(r, (j)). 
Then 

(1 d d 1 an ^ 

< r 1 \ = E^i 

2m [r dr dr J (3.5) 

$i(r, + 27r) = e-'^^^/'"'^^i{r, </>) . 

In this new gauge the electron wave function satisfies a free equation, but obey a 
non-trivial boundary condition upon making a trip around the solenoid. In other 
words, the Aharonov-Bohm effect is traded for the multi-valuedness of the wave 
function, the property anyons share. 

Indeed, the analogy is exact. Now we imagine that particles (in two spatial 
dimensions) have both charge e and flux /i, and that there exsit only charge-flux 
interactions, but no charge-charge or flux-flux interactions. Each particle, say a, 
creates a vector potential 

^fcj'x) = M - Xgj 

^ ' 27r|x-x„|2 

which is felt by other particles by the minimal coupling. Hence the Schrodinger 
equation is given by 



a=l 



$ = £; $ 

(3.6) 



_l4i=__£^V .7fc -Xb)k 
he " Stt^c^^ |x„-Xb|2 • 

This eqation is exactly the same as (2.6) upon identifying 

It is recognized that the anyon interaction is nothing but an Aharonov-Bohm 
effect. Since only the product of e and /u is relevant, one can phrase 

{charge : 1 
(3.8) 
flux : 2nce., 
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It's not exactly a Maxwell interaction, however, since there is no charge-charge 
interaction. We shall see in the next section that the Chern-Simons gauge theory 
precisely describes the anyon interaction. 

We would like to note that Aharonov-Bohm effects have wider applications. Let 
us consider the motion of a particle in an arbitrary multiply-connected manifold, 
A4. We can imagine that the space itself has nontrivial topology like x etc., or 
that the three-dimensional Eucledian space is obstructed by the presence of closed 
strings. Further we suppose that field strengths, or more specifically magnetic fields, 
Fjk = —djA'^ + dkA^, identically vanish in A4. Schrodinger equation is given by 

-^(\/ -i^aY ^0 = E ^0 inTW (3.9) 
2m V he / 

where djA'' — dkA^ = 0. In general Aj is not gauge equivalent to Aj = in a 
multiply-connected space. 

We define a new wave function by 

$i(x; C) = exp ( - ^ / dy A(y)l $o(x) (3.10) 

where C(xo,x) starts at xq and ends at x. For two paths, Ci(xo,x) and C2(xo,x), 
continuously deformable to each other, $i assumes the same value, since Fjk = 0; 

$i(x;(7i) =$i(x;C2) if Ci(xo,x) ~ C2(xo,x) . (3.11) 

Hence derivatives of $i with respect to x is well defined. 
$1 satisfies a free equation 

_^V2$i=£$i, (3.12) 
zm 

but is not single-valued. Let r(x) denote a transport of x along a closed path F. 
Then 

$i[F(x)] = w^(r) $i(x) 

(3.13) 



W^(F)=exp|-£^dyA(y)| 



If F can be continuosly shrunk to a point, then W = 1. W{T) depends on only 
homology of F with respect to M. It is a non-integrable phase factor, often called a 
Wilson line integral in the particle physics literature. We have seen that the gen- 
eral Aharonov-Bohm problem is traded for a free system with nontrivial boundary 
conditions. 

4. Chern-Simons Gauge Theory 

Consider a quantum field theory described by a Lagrangian 

^o=-^^'''''<^>^9-^P + i^^Do'^-^W? , ^^^^ 
Do = do+ iao , Dk = dk - ia!^ . 
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ci/n (flo = oP, ttk = —a'^) is a gauge field whose motion is described by the Chern- 
Simons term (cx e^^^'a^^yOp). (Rcfs. 63 70) V is a non-relativistic matter field, 
obeying either bose or fermi statistics. N determines the magnitude of the gauge 
coupling. A large |A''| corresponds to a weak coupling, as can be seen by rescaling 
Ufj,. We show that the system defined by (4.1) is equivalent to a neutral anyon fluid 
described by (2.6), and therefore by (2.4). 
Euler equations derived from (4.1) are 



N 
A-K 

idi 



1 (4-2) 

where fij,^ = d^a^ — d^a^ and 

i / \ (4-3) 

/ = -^(^^^feV'-(^feV')V) • 

The strengths /^^ of the Chern-Simons gauge fields are determined by the current. 
There is no physical degree of freedom for the gauge field. Hence the Chern-Simons 
field can be eliminated in favor of the matter field. 

It is convenient to take the radiation gauge div a = 0. The Chern-Simons field 
equation in (4.2) becomes 

— Aa*= = 

— Aao = dif - d2j^ . 

With an appropriate boundary condition, the equations can be solved to express 
in terms oi ■ 

In the case of a plane {B?) the boundary condition at infinity is subtle. A safe 
and rigorous derivation is obtained on a torus (T^). (Refs. 71- 85) We quote the 
result on a torus, taking the infinite volume limit. The argument presented here is 
very close to those in ref. 70 and in ref. 85. 

ao{x) = j dy hkix- y)j''{y) , 

(x) = {x) + (x) , 

, TTOe ik ^ -2 -1^ (4.5) 

ai(x) = -—— d'^Xk , ^{oia - d2a ) = ne 

iV ZTT 



(x) = j dy hj (x - y) {f [y) - rie) . 



Here rig is the average matter density, which generates a^{x). hj{x) is related to 
the two dimensional Green's function G(x) = (27r)~^ Inr by 

/.,(x)=-|e^'=a.G(x) = -^ . (4.6) 
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If a finite number of particles on an infintely large plane (i?^) are considered, 
one can set rig = in the above formulas. In applying to the superconductivity it 
is convenient to deal with a system with a finite density on R^. 

j^{y) in the expression for aP{x) in (4.5) contains a^{y) which is expressed in 
terms of if). Hence the Chern-Simons gauge fields are completely expressed in terms 
of the matter field. The substitution of (4.5) into (4.1) gives a Lagrangian which 
involves only ip and xp'^ : 

C\ = iip^ip — Hi 

1 , (4-7) 

Zm 

where a''{x) in D/. = dk — ia^ is given by (4.5). As it stands, the resultant Hamil- 
tonian Hi = J dx Hi involves four- and six-fermi interactions. There arises an 
ambiguity in ordering operators. It can be shown that the system defined by the 
Hamiltonian with the ordering adopted in (4.7) is equivalent to the system described 
by the Schrodinger equation (2.6). 

The Hamiltonian is not completely normal-ordered. It is given by 



Hi = ^ I dx [dki^^ + i^^a''){dk^ - iaV) 
2m J 

= Hi''>+Hi''>+H['^ 



(4.8) 



where 



- (DfcV)^(a;)(V'V(y) - ne)V(a;)} ^^'^^ 
Hi^ = ^ j dxdydz ft,fc(x - y)/ifc(x - z) 

and Dkip{x) = {d^ — id''{x))ip{x). Upon making use of J dy /ife(x — y) = 0, one can 
write /fp^ as 

= :T- f d^dydz hk{x-y)hk{x-z)i;^{x)i;Hy)i;Hz)iP{z)iP{y)i;{x) 
- j dxdy [/ifc(x - y)]VUa=)V'Uy)V'(y)V'(a;) • 



2m 



The equation derived from Hi is 



iij{x) = mx),Hi]^ICoi^{x) 

1 (4-11) 
/Co = --^^1 + oo(a;) + g{x) 
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where a''{x) in Dk and ao{x) arc given by (4.5) and 

^^"^^ = 2^ / ^^"^'^ ~ ^^^^^^ = 2^ J (^T^ ^^^^^^ • ^^-^^^ 

Eq. (4.11) differs from the classical Euler equation (4.2) by the g{x) term. The ad- 
ditional term is important to establish the equivalence between the anyon quantum 
mechanics and Chcrn-Simons gauge theory. 

Currents arc given by 



J°(x) = ijj^ijj{x) 

SHi SHi 

SVktp^x) ~ 5Vk^{x) 

-^{^^^ktp - {Dktp)^tP} (4.13) 



-^{i^\x)Dki^{x) - iDk^)\x)i;{x)} 



- — I dy hk{x - y) tp^x)il)^y)ip{y)il){x) . 



They are conserved: 9q J° + Vfe J*^ = 0. 

The Schrodinger wave function in quantum mechnics for a g'-particle system 
is related to the field operator ij) in (4.7) or (4.8) by 

$f(l,...,9) = (0|V.(l)---^A(rj)|vI/,) . (4.14) 

Here |0) and are vacuum and g-particle states, respectively, and tpifl) = ^{xa) 
where Xa = (i,Xa). For a system with a finite number of particles on a plane (i?^), 
one can put rie = and a'^(x) = in the above formulas. To obtain the Schrodinger 
equation, we differentiate with respect to t and make use of (4.11): 

»^ $f (1, ■■■,q)=f^ {0\m ■ ■ ■ li'ia) ■ ■ ■ V'((z)|*g) 

(4.15) 



^(0|V(l)---/CoV(a)---V'(9)l*,) 



a=l 



The following definitions and identities facilitate futher manipulations. First we 
define 

/Ci(a,6) = —hk{yia-^b){a''{xa)-a''{xb)} + '2g{xa,Xb) , 

TO 

V^{a,b) =^{Vl-Vl)hk{^a-^b) + ^hk{^a-^bf , 

TO TO (4-1d) 

1/3(0, 6, C) = —{hk{Xa-Xb)hk{Xa - ^c) 
TO 

+ hk{Xb - X„)/lfe(x6 - Xc) + hk{Xc - Xa)/lfe(Xc - Xfo)} . 
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Here 



9{x,y) = ^Jd^ hk{yi-z)hk{y - z) ip^^iz) . (4.17) 



/Ci(a, 6) is an operator, whereas V2{a,b) and V3(a, 6, c) are c- number functions. 
These operators with /Co(o) = K,o{xa) satisfy 



^(6)/Co(a) = {/Co(a) + /Ci(a, b) + ^2(0, 
V'(c)/Ci (a, 6) = {/Ci (a, 6) + 1/3 (a, 6, c)}V'(c) 



where 



^^^^ = E(V^)' + E EV.(a, 6) + ^ £ £ ^3(a, ^ c) 

0=1 a=2 6=1 a=3 6=2 c=l 



1 r 1 2 

z^^^Ei^^-^E^^^^--^")} • 



(4.18) 



and 

(0|/Ci(a,6)--- = . 
Apphcations of (4.18) and (4.19) lead to 

(0|^(l)---/CoVXa) ••>(<?) 1*9) 

a— 1 a— 1 b— 1 

= (0|{/Co(a) + J2 + ^2(«' + E E ^3('^' ^' ^)} 

6=1 6=2 c=i .4 20) 

X V(l)---V(9)|*,) 

= I ~ 2^^^"^^' + E ^2(a, 6) + £ y3(a, b, c) j $f . 

6=1 6=2 c=l 

Therefore the Schrodinger equation for is given by 

$f(l,---,g) = if(«) $f(l,---,g) (4.21) 



(4.22) 



a=l 67^0 

This is exactly Eq. (2.6), provided that 

Os = ^ (4.23) 

and that ipi^) satisfies anti-commutation relations. This establishes the equivalence 
of the anyon quantum mechanics and Chern-Simons gauge theory. The fermion 
representation is convenient to incorpolate the Pauli principle (2.3) for anyons. 
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5. Charged einyon fluid 

The anyon fluid described in the previous section is neutral. Anyons may be 
charged, interacting with each other electromagnetically. In the application to su- 
perconductivity, one needs to consider a charged anyon fluid. 

We have in mind material which has a layered structure as in newly discovered 
high Tc superconductors. The motion of electrons are mostly conflned in two- 
dimensional layers. The probability of the hopping of electrons from one layer to 
adjacent layers is very small. In many high Tc superconductors the resistivity of 
electrons in the direction perpendicular to CuO planes above Tc is 10^ to 10® times 
bigger than the in-plane resistivity. To the first approximation we may neglect the 
hopping interaction.*^ 

We shall adopt the "holon" picture originally advocated by Anderson. ^''"^^ In 
this picuture collective modes created by electron holes are spinless and charged. 
They are called holons, and are supposed to obey half-fermion statistics (^statistics = 
zt^n). The matter field denoted by iIj{x) corresponds to holon excitations. In 
our language ij^i^) satisfies anti-commutation relations, interacting through Chern- 
Simons gague fields with the coeflicient N = ±2 and through Maxwell fields with 
charge e. 

The electromagnetic interaction is not two-dimensional, however. Certainly 
there is a Coulomb interactions among electrons in two diflferent layers. Electro- 
magnetic waves can propagate in three dimensional space. 

In this article we consider two extreme limits. In one limit one can imagine an 
ultra-thin film which has a couple of, or just one, superconducting layers. Further 
we idealize the situation such that the only interaction among holons, other than 
the Chcrn-Simons or fractional statistics interaction, is the Coulomb interaction 
with a potential 1/r. We call it the ultra-thin film approximation. 

In the other limit we suppose material of an infinitely many layers (in the X1-X2 
plane) which are evenly separated with a distance d. Further we suppose that (1) 
electromagnetic fields Es, Bi, and B2 identically vanish, and (2) all fields Ei, E2, 
and S3 are uniform in the direction. Consequently all physical quantities such 
as the expectation values of currents (J'') are independent of 0:3. One can mimic 
it by considering a system in the idealized two-dimensional space, suppressing the 
third coordinate X3. It is called the two- dimensional approximation. Fluctuations 
of the El, E2, and fields are retained, in addition to the Coulomb interaction. 

We note that the two-dimensional approximation incorporates three-dimensional 
interactions in a specific way. Both the three-dimensional Coulomb interaction 
among electrons in distinct layers and the electron hopping between adjacent layers 
affect the three-dimensional motion of electrons or holons. When the electron or 
holon field is expanded in Fourier series in the xs direction with a momentum ks, the 
two-dimensional approximation amounts to retaining only the fcs = component. 

Real high Tc superconductors lie somewhere between the two approximation. It 
is necessary and important to have more thorough examinations of effects of the 
three-dimensional motion. No such analysis is available at the moment. 
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In the ultra-thin film approximation the Lagrangian is given by 

Do = do + i{ao + eAl^') , Dk = dk - iia!" + eA^^,) . 



(5.1) 



Here (/)(x) is an auxiliary field generating a 1/r Coulomb potential. A^^^ is an 
external electromagnetic field. £mm is bilinear in ip and ip'^ . After eliminating 
and (j), one obtains a Hamiltonian solely in terms of ij) and V'^- 



3D 

Coulomb 5 



(5.2) 



|x-y| 



V'(y)V'(a;) 



and flfc in Dk is given by (4.5). 

Equations of motion derived from Cfum are 

N 

= e{f - ne) (5.3) 

idoi> = [ - + ao + e{4> + Ag^^)} ^ 

where j'^ is given by (4.3) with the covariant derivatives in (5.1). 

In the two-dimensional approximation the Lagrangian is given, instead, by 

= --A-^e^''''a^^^ap + en,Aa + ^i^^Do^-^\Dki>\'' , (5.4) 

Do = do+ iiao + bAq) , Dk = dk - i{a'' + eA'') . 

Here the electromagnetic field contains both external and dynamical fields: A^i = 
A^^^ + A^y^. The corresponding Hamiltonian obtained by eliminating is 

H2D = J dx ^l-{Dki^)^Dk^)+eAoi>^i> + ^{El + B^)} (5.5) 

where Ek = -Fbfe and B = — F12. 

Equations of motion derived from (5.4) are 



N 

= ef - ene5^° (5.6) 
idoip = { - + (ao + e^o)} V' 
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Again j'' is given by (4.3) with Dk in (5.4). 

The Lagrangian forms Co in (4.1) (for neutral fluids) and £fiim / (for charged 
fluids) are convenient to develop a perturbation scheme. They are bihnear in ip or 
f/'^, and in the charged case the gauge invariance can be easily implemented in 
the perturbation scheme. On the other hand, the Hamiltonian forms Hi in (4.8) 
and -fffiim / H2B has an advantage that they involve only physical fields, and arc 
particularly suited to the evaluation of physical quantities beyond the perturbation 
scheme. We shall make use of both in subsequent sections. 

As we shall see, there is a subtle but important difference between neutral and 
charged anyon fluids. It seems that charged fluids (with a neutralizing background 
charge) are more stable than neutral fluids. 

6. Mean field ground state 

We consider an anyon system with a finite density rig 7^ on a plane. First 
we ask what would be the average Chern-Simons fields in the ground state. In the 
equation (4.2), (5.3) or (5.6), we replace the operator by its expectation value 
(j^) in the gromid state;. We expect that ()") = n,, and {j'^) = and all Maxwell 
fields vanish. In both neural and charged fluids we have 

6=^-^>^°^ (6.1) 

and /ofe = 0. In other words, particles, or holons in high superconductors, move 
in a uniform Chern-Simons magnetic field on the average. 

In the mean field approximation all gauge fields in (4.2), (5.3), or (5.6) are 
replaced by the average fields. The equation for the field operator tp is, in all cases, 

1 -2 
idotp = -—Dk V 
2m 

where the average vector potential is 

— e'^^ ^^^•^ = —e{N) e'^^ ^ in the symmetric gauge; 



a'' 



_ 27meX2 ^ _ ^ki ^ ti^g Landau gauge. 



(6.2) 



Here e{N) = +1 (—1) for A'' > (< 0), and the magnetic length, I, is defined by 

.0 lA^I 



27rne 

The corresponding one-particle Schrodinger equation 



(6.3) 



1 -2 

-^-Dfe Ua(x) = Ca Ua(x) (6.4) 
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is easily solved in either gauge. The energy spectrum is characterized by Landau 
levels: 

a = {n,p) n = 0, 1, 2, • • • 

(6.5) 

2)mP 

The first integer index n labels Landau levels. The second index p is cither a mo- 
mentum in the xi-direction in the Landau gauge, or an orbital angular momentum 
in the symmetric gauge. 
In the symmetric gauge 



ea= (n + 4)— TT = en 



s! 1 



,2 



(s + |p|)!27r/2 



(9 = 0,1,2,... and p e 2). 



(6.6) 



The Landau level index is n = s + 9{—p)\p\. Here L"{w) is the associated Laguerre 
polynomial 

i«(«;) = i «;-«e- £j(«;«+«e-) . (6.7) 

In the Landau gauge we impose a periodic boundary condition in the Xi- 
direction: Ua{xi + Li,X2) = Ua{xi,X2)- Then 



<^"d-"(x) = ^i= e-^'=-^^;„[(a;2-S2)//] , fc = ^ , X2 = e{N) kP 

(n = 0,1,2,... and pe 2). 
Here Vn{x) is related to the Hermite polynomial: 

^ ' (6.9) 



/CSO 
da; Vn{x)Vm{x) = 5n 



'U^p"^'*"(x) is a plane wave in the xi-direction, but is localized around X2 in the 
X2-direction. 

In a box (0 < X2 < L2), X2 must satisfy < 2:2 < ^2; or equivalently < p < 
LiL2/2ttI^. Hence the number of states per area for each Landau level, ul, is given 

by 

n.-^ . (6.10) 
Combinig (6.10) with (6.3), one finds that the filling factor, ly, is given by 

u=^ = \N\ . (6.11) 
riL 
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In other words, for an integer N, Landau levels are completely filled at least in the 
mean field approximation. It has to be stressed that this property holds irrespective 

of the density rig . 

If the tp field has spin i as electrons, the filling factor is given by v = ^\N\ so 
that the complete filling of Landau levels holds for an even integer N, provided that 
the magnetic moment interaction is sufficiently small as in the usual cases. We are 
mostly interested in the case |A^| = 2, corresponding to semions or half-fermions. 
There arises no qualitative change in physics properties, and we shall restrict 
ourselves in this article to spinless ip. 

We expand tjj in terms of {ua(x)}. 

V'(a;) = ^CaMa(x) , {Ca,c'/3} = 5al3 ■ (6.12) 
a 

The mean field ground state is given by, for = 2, 

I *mean> = H ^« ' 

aeG (6.13) 

G = {a = {n,p) ; n = 0,1} . 
The corresponding mean field energy is 

-Emean = mean I H['^ I *mean) = ^ • (6-14) 

aeG 

Here h[^^ is given in (4.9). Again putting the system in a box in the Landau gauge, 
one finds the energy density to be 

n=0,l p=l 

We recall that the energy density of a free spinless fermion fluid is exactly irnl/m. 
It is straightforward to check 

(j°(a;))^ean = (*mca„ I I ^-mean) 

= XI "a(x)'''Ua(x) 
aeG 

= Tie 

(/(^))mean = (*maan I {V^^fcV' ' (^fcV')V} I *mean) ^^'^^^ 

Y~ WaDkUa - {DkUa)^Ua} 



2m ^ 

aeG 



= 
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where = dk — ia''{x). 

7. Hartree-Fock ground state 

In this section we incorporate many-particle correlations in the Hartree-Fock 
(HF) approximation. The mean field approximation retains only h[^\ defined 
in (4.9), in the total Hamiltonian Hi, (4.8). In the Hartree-Fock approximation 

(2.) (3) 

"diagonal parts" of Hi and in (4.9) are taken into accoimt self-consistently, 
or equivalently to say, the ground state is determined to satisfy the Hartree-Fock 
equation. 

It was shown by Hanna, Laughlin, and Fetter^* that the Hartree-Fock ground 
state is exactly the same as the mean field ground state | 'I'mean), (6.13). 

I*Hf)= l*mean)= H ' " (^'l) 

However, it has a different energy. 

i^HF = (ifl)HF = (^HF I ffl I *Hf) ^ i^mean ■ (7-2) 

We first compute Shf- The original computation of Hanna, Laughlin, and Fetter 
was given in the first quantized theory. We present the evaluation in the second 

quantized theory. 

To facilitate the computations, we define the following sums. 



/(X,y) = ^ Ma(x)*Ma(y) =/(y,x)* 



aeG 



//s(x,y) = ^ Ua{xyiDkUa{y) = iDkf{x,y) 

aeG 



(7.3) 



We need to evaluate expectation values of various products of tl) and tjj^ . We 
denote (Q)hf = (*hf I Q I ^hf)- Then 

{^j;\x)^{x))uF = /(x,x) 

(^t(a;)^t(y)V,(y)V>(a;))HF = /(x,x)/(y,y) - /(x,y)/(y,x) 

(V'^ (x)^t (y)^t (z) V;(z)^(y)^(a;))Hp (7-4) 

= /(x, x)/(y, y)/(z, z) + /(x, y)/(y, z)/(z, x) + /(x, z)/(y, x)/(z, y) 

- /(x,x)/(y,z)/(z,y) - /(x, z)/(y, y)/(z, x) - /(x, y)/(y, x)/(z, z) 
Secondly 

i{'4}^{x)Dki}{x) - {Dktp{x))^tp{x))jiY = /fe(x,x) + /fc(x,x)* 
i{^\x)^^{y)i;{y)Dki^{x) - {D k^Pix))^ {yMy)i^{x))a^ (7-5) 
= {/fc(x,x) +/fe(x,x)*}/(y,y) - {/(x, y)/fc(y, x) + /(y, x)/fc(y, x)*} . 
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Making use of (4.9), (4.10), (7.4), and (7.5), one finds 

£uF = -^{{h['^ + fff ' + i^f ^)>HF = + + (7-6) 

vol 

where vol is the volume. £^^^ is the same as the mean field energy density 

cil) _ c ^ !!Zi /7 7^ 

m 

and 

= '^'^'^y /ife(x- y)[ {/fc(x,x) + /fe(x,x)*}{/(y,y) - rie} 

- {/(x, y)/fe(y, x) + /(y, x)/fe(y, x)*} " 
f = — -— / dxdydz hk{x - y)/ife(x - z) 

2/(x,y)/(y,z)/(z,x) - /(x,x)/(y,z)/(z,y) 

- 2|/(x, y)|2{/(z, z) - n,} + /(x, x){/(y, y) - ne}{/(z, z) - rie} 
+ / '^'''^y ['^fe('^-y)]' {/(x,x)/(y,y)- |/(x,y)p} 

The quantities /(x, y) and /^(x, y) in (7.3) depend on the gauge chosen. The 
symmetric and Landau gauge defined in (6.2) are related by 

asym(x) = aLndau " VfcA(x) 

X1X2 (7.9) 



(7i 



A(x)=e(iV) 



2/2 



If ■u|^™(x) is a sohition to the Schrodinger equation in the symmetric gauge, then 
g«A(x) (x) = (x) is a solution in the Landau gauge with the same energy 

eigenvalue. It is a linear combination of M^p"^^"(x) in (6.8). With a given Landau 
level (energy eigenvalue), the sets {^^^"'^'^'^(x)} and {u^p'"^^"(x)} are related by a 
unitary transformation. 

Hence, if the set G in (6.13) represents completely filled Landau levels as in the 
case under consideration, then 

^ t.^andau^^)*^Landau^y^ ^ ^ ^Landau ^Landau 

so that 

/(x, y)^'^- = e^^^(-)-^('')> /(x, y)L-d- 
/fc(x,y)^'''" = e^^^W-^^'^)} /fc(x,y)L-dau _ 
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Thanks to the relation (7.10), every term in (7.8) is separately gauge independent. 

It is easiest to evaluate /(x, y) in the Landau gauge. Making use of (6.8), one 
finds 

\N\-1 oo 

/(x,y)L-dau = 5^ ^ — e^'=(--^^) V4{X2 - X2)/l] Vn[{y2 ' y2)/l] 
n=0 p=-oo 1 
. I^l-l .oo 

= ^ E / e^'^-^-y^^ Vn[{X2 - X2)/l] Vn[{y2 - ^2)/^] • 

Employing explicit forms of u„'s and integrating over fc, one obtains 
TV = ±1 : /(X, y)— = ne • exp { - 1^ ± ^ ^^^ - V^)^-^ + V^) 

= /o^(x,y) (7.12) 
N = ±2 : /(x,y)"- = |l- (^il^l./±(x,y) 



Note that l"^ = \N\/2'Kne. fki^,y) is obtained from (7.3). 

N = ±l : /i(x,y)L-<iau^_L|(^^_y^)^i(^2-y2)}/o^(x,y) 

/2(x,y)L-<iau ^ ^•y^(^^y)Landau 



7V = ±2 : /fe(x,y)L^"'i^'^ = |l - \ ■ h{^,y)'%-^f^ 



(7.13) 



21^ 

We are ready to evaluate (7.8). We drop the superscripts 'Landau' in /(x, y) 
and /fc(x, y). Note that /(x, x) = rie and /fe(x, x) = 0. Furthermore, 

/ife(x-y)/fc(x,y) /(x,y) . (7.14) 



Therefore S^"^^ becomes 



^ for7V = ±l 



for AT = ±2. 
4m 



For 5(3) 

we have 

1 1 f ^ ^ ^ 1 (x-y)(x-z) 



vol 2m J 



N2 (x-y)2(x-z)2 
'2/(x,y)/(y,z)/(z,x)-ne|/(y,z)|2 



(7.16) 
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As will been seen below, £2^^ and £"3^^ have infra-red divergences, which cancell each 
other. 

In evaluating we rename y — x and z — x to be new y and z, respectively. 



Then 



^^^^ " J ^^'^^ fi'Jv(y)3Jv(z)ffjv(y - z 



y z 

^ '^''P { " + z^ + (y - z)2] ± ^(yi^2 - y22i)} ^j yj-^ 



r 1 for TV = ±1 

ff^W= foriV = ±2 

We introduce polar coordinates by 

The integrand is independent of ^, and yz/2Z^ = rpcoscj) and — y2Zi)/2P 
rpsin^. (7.17) is transformed to 

3 poo p27T 

^ = —f- ■ 27r • 2/2 / drdp / d^i. cos ^ • exp | - - + ^pe^^^"^ | 
^^rn Jo Jo 

1 for AT = ±1 

lr2)(l - ip2)(l - ir^ - i 

Integrating over (f), one finds 
N^m Jn 



(1 - ir2)(l - ip2)(i _ 1^2 _ 1^2 ^ ^pcos0) for A = ±2 



TTrp 



7rrp(l - ir2)(l - lp^)[2 - l{r^ + p^) + p^^ 



^ forA = ±l 

2m 

^ for A = ±2 

8m 



(7.18) 



Similarly, the second term, £3^^ in (7.16) becomes 

= -t4^ / dyd. 4^ ,.(y - z)2 e-(-^)V-^ 
2mN^ J y^ z^ 

/■ , , y(y-u) . x2 -uV2i^ 



(7.19) 



2mA2 

where we have introduced u = y — z. This time we define 

yi+iy2 = V2l re'^ , ui+iu2 = V2l pe'^^+'^^ . 
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Then 



'(3) 



277/ 



2„3 



The ^-integral gives 



p r 
drdpdcj) - 



rp cos ( 



J, j,2 _|_ p2 _ 2rp cos ( 



QNinf e-" 



rp cos ( 



2Tre{r - p) 



10 + p2 — 2rp cos 

The rest of the computation is straightforward. We find 



:-(3) 



mm 



dr — 



+ 



The integral in (7.21) diverges in the upper limit. 



The evaluation of fg^' is easy. 



ci3) 

^3 



n: 



2TOiV2 

mm 



dr 



N,±2 



The divergent integrals in (7.21) and (7.22) cancell each other. 
Adding (7.18), (7.21), and (7.22), one finds 



1 nni 



2 m 
5 im'i 



for N = ±1 
for N = ±2 



f(3) = 

32 m 

Finally combining (7.7), (7.15), and (7.23), we obtain 

1 7rn3 



-HF 



V 2/ m 2m 

/ _ 1 ^\ 7rn| 29 Trn^ 



32 m 



for iV = ±1 
for = ±2 



(7.20) 



(7.21) 



(7.22) 



(7.23) 



(7.24) 



The correction to the mean field energy is large for A = ±1, but is relatively small 
for N = ±2. 

There are various ways of showing that | ^'hf) in (7.1) is the Hartrcc-Fock 
ground state. The Hartree-Fock approximation amounts to finding an approximate 
ground state for a many-particle system in the form of a Slater determinant formed 
from one-particle wave functions. In the language of the second quantized theory 
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we expand the field operator ip{x) in a complete orthonormal set {mq(x)} to be yet 
determined, and write a trial ground state as 



y dx {ij)^ ll){x)) = He ■ vol 



where (Q) = (\I>g | Q \ '^g)- {^ai^)} and the set G are determined to minimize the 
expectation value of the Hamiltonian {H). 
In Eq. (4.11) 



ii'ix) = I - + ao{x) + g{x)^ i,{x) (4.11) 

the approximation amounts to retaining only diagonal pieces on the right hand side. 
For instance 

'4)\y)il)\z)^{z)^{y)i){x) (7.26) 
i^PHyWy)) {4'Hz)i;{z)) ^Pix) - {4>\yMz)) (V^(z)V(j/)) V'(a:) 

- (ipHyW^)) {^Hzmz)) i^{y) + {ipHyMz)) {^Hz)i;{x)) ^{y) 

- {^\ymy)) {^\z)^{x)) ^{z) + {^\y)^{x)) {^\z)i^{y)) ^{z) . 

The set {mq(x) ; a E G} is determined to satisfy Eq. (4.11) with the above substi- 
tution made. The equation thus obtained is called the Hartree-Fock equation. 

It was shown by Hanna, Laughlin, and Fetter^^ that | completely filled Landau 
levels formed from the mean field eigenstates (6.6) or (6.8) satisfy the Hartree-Fock 
equation, and therefore | 4'hf) in (7.1) is the Hartrcc-Fock ground state. The 
computation is similar to, but more complicated than, that of (ffi)jjp presented 
above. Readers should refer to the original paper for details. 

The Hartree-Fock equation can be written solely in terms of the fimction /(x, y) 
defined in (7.3). With the ansatz (7.25), the expectation value {H) can be viewed 
as a function of (^/^'''(x)V'(y)) = /(x, y). 

(iIi)=H[/(x,y)] , /(x,y)*=/(y,x) . (7.27) 

/(x, y) is determined by the extremum condition: 

S'M. f 
-^j^ y = subject to J dx/(x, x) = rig • vol . (7.28) 
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/(x, y) is determined by the equation up to an arbitrariness due to gauge degrees 
of freedom. Eq. (7.28) is equivalent to the Hartree-Fock equation. 

8. RPA and SCF 

There are approximation schemes which he between the mean filed and Hartree- 
Fock approximations, and are useful to investigate various physical quantities such 
as the excitation spectrum, current-current correlation functions, response to exter- 
nal perturbations, and so on. They are the random phase approximation (RPA), 
self-consistent field method (SCF), and hydrodynamic description. In the anyon 
model under consideration, all these three are equivalent to each other. 

In a neutral anyon fluid RPA is formulated from the Hamiltonian Hi[iIj,iIj'^] 
in (4.8), whereas SCF from the Lagrangian jCo[tp,tp^ jU/j] in (4.1). The difference 
lies in whether Chern-Simons fields are integrated out first, or retained. A similar 
statement is valid for a charged anyon fluid. 

For definiteness we shall restrict ourselves to a neutral anyon fluid in this section. 
One can establish the diagram method, or Feynman rules, from the Hamiltonian 
Hi. ff^^^ in (4.9) defines a bare propagator for the ^-field. It is a propagator in 
the mean field, depicted by a solid arrowed line. and define interaction 
vertices. 

(2) 

The two-body interaction generated by Hi is given by Va in Fig. 3. A dashed 
line corresponds to /ifc(x — y), representing virtual "propagation" of Chern-Simons 

fields. A crossed circle at the vertex x indicates the derivative factor iDk- 

For the purpose of establishing Feynman rules, it is convenient to start with the 
normal-ordered form (4.10) for h['^K It yields three- and two-body interactions VJ, 
and Vc in Fig. 3, respectively. Both involve two dashed lines. 



Fig. 3 Feynman rules derived from (4.8). The rules are simplified by re- 
taining Chern-Simons fields as independent variables. See Fig. 7 in 

Section 10. 

There are important rules resulting from the form of the Hamiltonian. Recall 

that {'iIj'^ {x)tp{x)) = He, {tp^ {x)D ki^{x)) = 0, and / (ix/ife(x) = 0. Therefore, con- 
traction of two solid lines at the same vertex in Va yield a vanishing result. Similarly, 
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contraction of two solid lines at either y or z vertex in Vb yields a vanishing result. 
See Fig. 4a. However, contraction of two solid lines at the x vertex in V;,, and at 
the X and y vertices in Vc yield a non-vanishing result. See Fig. 4b. 



Fig. 4 Feynman rules — constraints. 

Let's consider correlation functions of the currents J^^{x) defined in (4.13). Note 
that J''{x) defines two- and four-point vertices. (See Fig. 5.) To all order in 
perturbation theory 

( J^(:r)) = ne (5'^° . (8.1) 



Fig. 5 Currents J'^{x) . 

For the correlation function for a neutral fluid 

D!:i%x,y) = -i{T[P{x)J^{y)]) , J^{x) = J'^ix) - 5^° (8.2) 

RPA constitutes in keeping only daigrams in which no dashed line is a part of closed 
loops involved. Typical diagrams involved are depicted in Fig. 6. In RPA we 
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have only sums of bubble diagrams. Further, the four-point vertex in J*^(a;) does 
not enter at this level. 



Fig. 6 r>5J^ 



'{x,y) in RPA. 



One can write down Schwinger-Dyson equations for D^^^ , which was done by 
Fetter, Hanna, and Laughlin,^^ and by Chen, Wilczek, Witten, and Halperin.^^ Due 
to the presence of the off-diagonal couplings, the equation takes the form of a 3 x 3 
matrix. In the following section we shall show an alternative method to evaluate 



The self- consistent field method starts with the field equation (4.2), iretaining 
the Chern-Simons gauge fields:^^'^^ 



We suppose first that there is a consistent gauge field configuration a^{x) which 
is treated classically. Secondly, with this given ap(x) we determine a quantum 
mechanical state vector, | \E'[a^]), for the matter field tp{x), by solving the second 
equation of (4.2). Thirdly, we replace the current on the right side of the first 
equation of (4.2) by its expectation value: 



477 ^ 



(4.2) 




(8.3) 
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Eq. (8.3) is solved for a^(a;) to find a self-consistent field configuration. 

We have known that the configuration, (6.2), of the mean field ground state, 
which is the same as the Hartrcc-Fock ground state, solves Eq. (8.3). We arc socking 
for more general, x-dependent solutions. If a deviation of a^(a;) from (6.2) is small, 
one can employ a perturbation theory to find a solution. 

For timc-indcpcndcnt configurations the procedure is particularly simple, as 
was first performed by Hosotani and Chakravarty.'^^ The generalization to finite 
temperature is done by Randjbar-Daemi, Salam, and Strathdee,^* and by Hetrick, 
Hosotani, and Lec.^^ Time-dependent configurations have been recently analysed 
by this method by Chakravarty.^* 

Consider a time-independent, small fluctuation of Chern-Simons gauge fields, 
a'j}\x) = an{x) — ai^{x). We first solve the one-particle Schrodinger equation with 
this an{x): 

{-^^fc + «o(a;)| ;/a(x;aW) = e„(aW)u„(x;a«) . (8.4) 

Both ea(a^^^) and Ua{x;a'j}^) are determined perturbatively. {Ua{x;a^^^)} defines 
a complete, orthonormal basis, with which we expand tlj{x) as 

^(a;) =5^c„(aW) w„(x;o«) , {c„,4} = ^„/3 . (8.5) 

a 

So long as alj^"* (x) is sufficiently small and smooth, the spectrum retains the 
structure of Landau levels, although they are not degenerate any more in general. 
With the same set as G in (6.13), we define a state 

vl>[aW]= n4(a^'^)|0) , (8.6) 

from which the current in (8.3) is determined as 

J^[x;a^]=neS>''>- j d'yr^''{x,y)ai'\y) + --- . (8.7) 

At finite temperature T(= /3~^) we evaluate the matter part of the free energy 
with a given a^(a;) by 

g-/3-PKl = Xr e"''-^"!''''"^'! (8.8) 



where 

Ho[^,af,] = j 1^ (i?feV)^(^feV') +«oV'V 



(8.9) 



The current is given by 

Sa^{x) 



Sa„{x) 
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which leads to an expression similar to (8.7). Incorporation of electromagnetic 
interactions is straightforward. 

In this section wc have explained the two approximation methods, RPA and 
SCF. These two look quite different from each other. RPA is defined for Green's 
functions for matter fields in terms of Feynman diagrams, whereas SCF is written 
in the form of gauge field equations. Wc shall show in the next section that these 
two are indeed the same, and are equivalent. 

9. Path integral representation 

A bridge between RPA and SCF becomes most transparent in the path integral 
formalism. One can deal with both neutral and charged anyon fluids at once. The 
key step is to consider the transition amplitude or partition function in the presence 
of external gauge potentials. 

To simplify notations, we write as 



(9.1) 



where Dq = do + iag and Dk = dk — ia^ ■ In terms of these definitions 

£2D a. A] = [a] + C^^^ [A\+Cf[t,a + eA] 

Consider a charged anyon fluid described by (5.4) or (5.5). Wc introduce an 
external gauge potential a°^'(x) in a gauge invariant way, making a replacement 

Cf[i);a + eA]^ Cf[i);a + eA + a''''^] . (9.2) 

In other words, the external Lagrangian is given by 

jCe^t= Cf[ip;a + eA + a^''*]- Cf[i:;a + eA] 

— "o J + "extJ 2j7j v"extJ J 



(9.4) 



where 

f{x) = ^p''^ip 

3\x) = -^(V^tVfeV - VfeV^ • V) - ^V'V • (a'^ + eA'=) 
In the presence of an external potential a total gauge-invariant current is 

^ ^ (9-5) 

itot=i° > itot = / - - OextV-V • 
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The transition amplitude at T = is given, in tlie patli integral representation, 

by 



exp(a2DK''*]) = JvAVaViP^Vi; expji J dx (r^*^ [a] + [a] 
+ C^^ [A] + Clf [A]+Cf[t,a + eA + a^^'] ) | 



(9.6) 



Here £g f [a] and are the gauge-fixing terms for and A^, respectively, 

and shall be specified shortly. The initial and final state wave functions have been 
absorbed in the definition of the path integration measure. Their explicit forms are 
irrelevant to compute various quantities described below. 
Several identities follow from (9.6). First 

"5a-Hx) = *5a-t(x) "^''^ = • ^^-"^^ 

Recalling that {j^^t )oext=o ~ ^gS^^, one finds a current dynamically induced to be 







SI.. ,0 . (9-8) 



Sa^^ {x) 

Similarly the second derivative leads to 
(5a^'^*(x)(5ay''*(t/) 

= iT\jUx)jUy)])a... + i(l - S'^')5^^ ^ ifix)),^^^ S{x - y) (9.9) 

<5^-?2D _|_ <5/2D SI2D 



In the a^'^*(a;) = limit one has 

{Ty^{x)f{y)]),^^^^o-S^°6''° nl 



6a^^*{x)Sal^^{y) 



- i(l - (5''°)(5^" — S{x - y) . 

aext=0 ™ 



(9.10) 



Formulas for a neutral anyon fluid arc obtained from the expressions above, 
dropping electromagnetic fields A^ entirely. The amplitude becomes 



exp(i/nK'^*]) = j VaVip^V^ 

X exp |z J dx (£^^[a] + C';i[a] +Cf[t,a + 0*"=*]) | 



(9.11) 
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The induced current and correlation function are given by 

{Tmx)f{ym^'=:^ - 5^'5^' nl (9.12) 



i(l - ^"0)5'^'' ^ b[x - y) . 



aoxt=0 



The first equation of (9.12) is related to (8.7) in SCF, whereas the second equa- 
tion gives the correlation function (8.2) in RPA. Both are derived from the effective 
action J„[o°'^*] or /2D[a'''^*] for a neutral or charged anyon fluid, respectively. The 
connection between RPA and SCF is established by evaluating the effective action. 

As for gauge-fixing, it is most convenient to take 

(9.13) 

Cft\A] = -{d,A^f 

gives the standard Lorentz covariant Feynman gauge for electromagnetic 
fields. We have retained the gauge parameter a for Chcrn-Simons fields. For a = 1, 
£g f [ a ] gives a spatial Feynman gauge. In the a = limit it reproduces the 
radiation gauge Vfea*^ = 0. 

For Chern-Simons gauge fields we have 



= neaW + laW (Aq + Ag.f,)''- 



(9.14) 



where the kernels Aq and Ag.f. arc given by 

, Ag.f. =-- I dl did2 I . (9.15) 





Note that a{,^^ = (4^\o^^\4^^) = (a(i)°, -o(i)\ -a^^^^^. Integration by parts has 
been made in the above formulas. The Chern-Simons term itself is singular, as 
detAo = 0. However the total kernel is regular: 



A = Ao + Ag.f. 
For electromagnetic fields we have 



(9.16) 



1 



£™ [A] + /:™ [A] = -A^ g'^'^d^ A, + en.A^ (9.17) 
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where the metric is g^"^ = diag (1,-1,-1). 

The integration over various fields in the formula (9.6) can be done in any order. 
Integrating first over the Chern-Simons fields a^(.x) in the limit a = is equivalent 
to eliminating them to get the Hamiltonian (4.8) or (5.5). An alternative is to 
integrate the matter fields ip{x) and ip^x) first, maintaining the symmetry of the 
CS and EM gauge coupHngs. 

To see that the Hamiltonian (4.8) is reproduced in the neutral case by integrating 
Ufi, we write 

1 aWA'*''a(i) + Lfli,; a] = ^ a«(A + Era^J^ + a^J'' + a^] (9.18) 
where 

S^'^ = -(l-,5^°)(5''^— 
m 

? =j°=VV (9.19) 



Further we note that 



A-' = -— -c% I + 0(a) 



+^1 / (9.20) 



N V2 

(A + S)-i = A-i + A-^SA-i + A-^SA-^SA-^ + ■ 
In the a ^ limit 



A"^SA"^E:A"M„=o = etc. 




so that the integration over a|/^ yields a Lagrangian 

= J-f'^YjO ^ -0 ^ aO kl --k 

2m\N ) ' ' ■' N V2 



Noticing that 



/ 



(9.21) 



we observe that (9.21) yields the same Hamiltonian as in (4.8). The path integral 
formalism is equivalent to the operator formalism. 
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10. RPA = linearized SCF 

In practice it is more convenient to integrate the fermion fields ip and ip^ first 
to evaluate 7„[a**''*] in (9.12) or 72D[a'*''*] in (9.6). Let us define 

exp{iSf[a\} = J Vip^Vil) expi^i J dxjCf['ip;a\^ . (10.1) 

It gives the transition amplitude for the fermion fields in the presence of gauge fields 
a^(a;). Then we have 

exp{i/n[a^''*]} = J Va exp l^iSf[a + a"""^] + i J dx {jC^^[a] + jC^i[a])^ 
exp{z/2DK''*]} = J VA expl^il^eA + a'"''] (10.2) 

+ i I dx{Cf''[A]+C^f,[A])'^ 

To evaluate Sf[a] we decompose C flip; a] into the zeroth order and interaction 
parts: 

£/[V';a] =£/[V';a(°)]+4"* . (10.3) 
Here a{j°^ (x) is the average field configuration given in (6.2) and 



^Wti)fcV'-(^feV')V} (10.4) 

_(a(i)'=)2^^t^ 
^ ^ 2m ^ ^ 

where a^}^ = - a^°^ 

Cf[ip; a^°^ defines a propagator of the ip field in the background potential a^"', 
whereas defines interaction vertices containing a'l^\x). The first and second 
terms in (10.4) give one gauge field leg, whereas the last term gives two legs, as 
depicted in Fig. 7. 

Sf[a\ is nothing but the ofFcctivc action for a,fj_{x) generated by dynamics of ijj 
and f/'^ fields. The standard diagram technique can be employed. Since tp and 
are integrated, fermion lines must be closed. Sf[a] — Sfla'^^^] is the sum of connected 
diagrams. Further, since the interaction £j"* is bilinear in tp and tp\ diagrams thus 
generated are all one-loop. One can arrange them according to the number of legs 
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of gauge field o^^^ as in Fig. 8. 



Fig. 7. Vertices generated by in (10.4). "0" and "k" at the ends of 

(1) (l)k 
dashed lines indicate Oq and w- ' , respectively. 



Fig. 8. The eflfective action Sf[a]. 

Contributions coming from diagrams (a) and (b) in Fig. 8 are easy to evaluate. 

diagram (a) = —a^^^ (V^V-') = — ^o^^ rie 

1 (10.5) 
diagram (b) = - — (a'^^^*')^ne . 

Note that in the nonrelativistic system under consideration we always have 
{ip^'ip{x)), instead of limx^y{—l){T[il){x)il)^{y)]), in the first-order perturbation. A 
similar diagram containing the second vertex in (10.4) vanishes, since {j'^) = 0. 



32 



Neutral and Charged Anyon Fluids 



In general one has an expansion 

Sf[a\ = Sf[a^"'^] - I dxneai\x) 

(10.6) 

+ j dxdy - a«(a;) T'^''{x - y) a<i}Hy) + ■■■ . 

In writing the T'^'^ term, we have employed the translation invariance of the zeroth- 
order system described by Cf[il),a^^^]. In the momentum space it becomes 

/ ^ ^<^H--,-q)r-(u.,q)«W(c.,q) . (10.7) 

Contributions from diagrams (c), (d), and (e) to F'''' will be evaluated in the fol- 
lowing sections. 

The next step is to integrate over Chern-Simons fields a^(a;) in (10.2). Again a 
diagram method may be developed for the integral. Higher order terms in (10.6), 
namely terms involving three or more a^^'^'s, give "interaction" vertices. 

By droping all these higher order terms the system is linearlized. We call it "ifte 
linear approximation^^ . The resulting integral is a simple Gaussian integral, whose 
evaluation is straightforward. As we shall see shortly, RPA and the linearlized SCF 
are nothing but the linear approximation. 

Let us consider a neutral anyon fluid. We observe 

Sf [a + a^^^*] + Jdx (£^s [^j + £cs [^j) 

= Sf[a^'n -Jne (4^) + ) + J i(aW + a-*)r'''^(a(i) + a^) + ■■■ 
+ (|nea«+/ia«A-aW) 

+ J l{a^^^ + a'"* r (A + r)-i} (A + T) {a^^) + (A + T)-^T a^'^'} 
- j ia'^''* r (A + r)-ir a'^'^* + j ia"''^* T a'^'^* . 

(10.8) 

We have suppressed a measure dx or dudq in the experssion. The integration over 
a'jj^ immediately leads to 



^ - I dxn, a'o'^'ix) - / dxdy -a^^^^) Q^'i^ " v) aTiv) + ' ' ' (10.9) 



where 

Q„ = r(A + r)-ir-r 



1 1 

(10.10) 



i + A-ir l + TA-i 
= - ^ + ^A-l^-^A-l^A-l^ + • 
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represents a propagator for a{i^^ Therefore Q„ has a diagram representation 
given in Fig. 9. 



Fig. 9. Diagrams for Qn- Diagrams for Qc are obtained by replacing each Chern-Simons 
field propagator (dashed line) by the sum of propagators of Chern-Simons fields 
and electromagnetic fields. Compare (10.10) and (10.20). 



From (9.12) the jj-correlation function is given by 



(10.11) 

- Q'n'ix - y) - (1 - S>^'>)6>^'' ::is{x-y) + --- . 

m 

We recall that the contribution of diagram (b) in Fig. 8 to F is given by (10.5) so 
that 

^T%^=P'^S{x-y) , (10.12) 

TO 

which is precisely the negative of the last term in (10.11). In a normal metal this 
is the end of cancellation. In anyon fluids something special happens. As shall be 
shown in section 12, diagram (c) in Fig. 8, in part, yields the same contribution as 
(10.12) with the opposite sign. Hence the last term in (10.11) survives. 

The series generated in (10.11) with (10.10) substituted exactly reproduces the 
diagrams in Fig. 6 in Section 8. Hence 



= (x - y) - (1 - SnS''-' - 6{x - y) . 



,o^x^. »e . (10.13) 

TO 

RPA is equivalent to droping, in 7n[a''''*], all terms cubic or higher-order in a^''*(a;). 

34 



Neutral and Charged Anyon Fluids 



How about SCF? We return to the expression (8.7). In SCF one determines 

the current in the presence of non-trivial gauge fields 0^(2;), which is nothing but 
evaluating Sf[a] in (10.1). Since j''(x) = —6 J a]/(5a^(x), one immediately 
finds 

6Sf[a] 



(10.14) 



We recognize that T^^" appearing in (8.7) is the same as F'^" defined in (10.6). 

In SCF the current thus obtained is inserted into the Chern-Simons field equa- 
tion (4.2), which is then solved to determine a sclf-consistcnt nontrivial field config- 
uration. If higher-order terms in (10.14) arc dropped, the resulting field equation 
becomes linear. Solving the equation is equivalent to performing a Gaussian inte- 
gral in the path integral formalism, since the latter amounts to picking a stationary 
path of the action. 

An alternative way of seeing this is to examine a response to an external field 
in SCF. Since Jext = Aoext, the field equation in SCF, with the gauge fixing term 
added, becomes 

A a^^^ = Jind + Jext = (-r a(i) + ■■■)+ Attest , 
from which it follows that 

j.^yinearizedSCF ^ _r (A + Tr'Aa,^^ = -F (1 + A'^r)"! flext 

= Qnacxt (10.15) 

T linear 

— ^ind 

We have stressed in the last equality that the expression is exactly what one obtains 
from the first equation of (9.12) combined with (10.9). The relations (10.13) and 
(10.15) together establish the equivalence between RPA and the linearized SCF. 

The generalization to a charged anyon fluid is easy. In the two-dimensional 
approximation one needs to do one more integration over A^{x) in (10.2). With 
the expressions (10.9) and (9.17) inserted, the exponent of the integrand becomes 

-HeieAo + ag"') - ^ (e^ + a^"*) Qn {eA + a^^') + ■■■+ ^APA + eUeAo 

= 1 {A - a-* eQn {P - e^Qn)"'} (P - e^Qn){A - (P - e^Q^)"^ eQ„ a^*} 

- I a-'"* eQn {P - e^Qn)-' eQn a''^' - i a^^' Qn a''^' - riea'^' + • • • 

(10.16) 

where 

pi^^ = gl^'^ Q'^ (10.17) 
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Therefore we have 



/2dK"*] = -jdx neal^\x) - j dxdy i a-*(x)Qr - y)aT\y) + ■■■ (10.18) 
where 

^ l-Qn e^P-i ^" = ^" 1 - e2p-i Q„ (10-19) 

— ~1~ Qn ^ -P ~1~ Qn ^ Qn^ ^ Qn ~1~ * * * 

Combining (10.10) and (10.19), one finds 

„ 1 1 „ 



Qc 



1 + (A-i + e2p-i)r l + r(A-i+e2p-i) 



(10.20) 



The proper vertex — F, which summarizes one-loop fermion interactions, is con- 
nected to the next one by a propagator of either Chern-Simons fields (A~^) or 
electromagnetic fields (e^P^^). There is only one proper vertex, since both Chcrn- 
Simons and electromagnetic fields minimally couple to the fermions. The final 
expression for Qc above is obvious from the viewpoint of the diagram method. (See 
Fig. 9.) 

To summarize, in the linear approximation, which is equivalent to RPA and the 
linearized SCF, 

(x, y)^'"^^^ = Qr - y) - (1 - ^''°)5''^ ^ - v) 

T linear 



11. Response function 

As it stands from (10.13), (10.15), and (10.21), Qn or Qc determines a linear 
response to an external perturbation. It defines a response function. Location 
of poles in the response function Q„ or Qc gives an energy spectrum of particle- 
hole excitations. In passing, location of poles in the fermion propagator yields a 
spectrum for fermionic excitations, which is quite different from that of particle- hole 
excitations in anyon fluids under consideration. 

One can also examine a response to an external magnetic field in a charged anyon 
fluid, from which the cxsitcncc or non-existence of a Mcissncr effect is checked. 
Conductivity or resistivity tensors can be computed from Qn or Qc- Examinig a 
response to external density perturbation Jext(^) gives information on muon spin 
relaxation. Many other consequences can be drawn from Qn or Qc- 
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We have seen in the previous section that Q„ and Qc are determined by the 
proper vertex F defined in (10.6). There are many restrictions resulting from the 
definition (10.6) and the conservation law. 

First, from the definition wc have F^'^(x — y) — T'^^{y — x), or 

F''''(a;,q)=r'''^(-a;,-q) . (11.1) 

Secondly, the gauge invariance or current conservation implies d^T^^'^^x — y) = 
and dyTt'^ix -y)=0, or 

q^T^-'{uj,cD=0 = q.r^''{uj,ci) (11.2) 
where = (w, — q). Thirdly, the rotational invariance implies that 

roo = A 

=qjB + e,kqkC 
=qjB' + e,kqkC' 
P'' = 6jkD + ejkE + qjqkF 

where A ^ F are functions of uj and q = |q|. 

Relations (11.1), (11-2), and (11.3) lead to a decomposition^'^'^^ 

F°°(c.,q)=g2no 
F°-' (w, q) = ujqjUo - ie^kqkT^i 
r-^°(a;,q) = wq^IIo +iejkqkT^i 

r^''{u>, q) = 5jk w^rio + icjkUiUi - {q'^Sjk - qjqk) ^2 



(11.3) 



(11.4) 



where all II^'s are functions of and q^ only. If the perturbative ground state is 
stable, Sf[a] defined in (10.1) is real, and therefore from (10.6) r'''^(a;)* = r'''^(a;), 
or T'^'i-uj. -q)* = F^''(u;,q). In other words, Hfe's {k = 0, 1,2) are real: H*^ = H^. 
In a frame q = {q, 0), 



/ g^FIo uqUo 




1 1 w^Ho +iujUi I . (11.5) 

To evaluate Q„, we recall Eq. (10.10): 

Q„ = r(A + F)-iF-F 

= -A(A + F)-iF (11.6) 
= -F(A + F)-iA . 

Take a frame in which q = {q, 0). F is given by (11.5), and 



(11.7) 











io 


1 


• 
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Combinig (11.5) and (11-7), one finds 



A + r = 



ujqllo +^9111 



\— iglli —iwYii 112 / 



{ll.i 



- ^ N 

Hi = Hi , 

^ ^ 27r ' 

112 = w^Ho - g^Ha . 
We need to evaluate (A + r)~^. First we note 

fl* - -2 

det (A + r) = ^(noHa-nr) 

A straightforward manipulation leads to 

(A + r)-i = 
1 



(11.9) 



/ Ha 



g2(non2-ni) 



2 -iglli 





(11.10) 



Both A and (A + r)-^ depend on a, but the product A (A + T)"^ or (A + T)-'^K 
is independent of a. (We shall see shortly that the current conservation guarantees 
the a-independence.) Q„ is given by 



non2 - Hi 



















fia / 


N 


/ 



-iq\ 
—ioj 




V +iq +iui 


/ 



(11.11) 



This expression was first given by Aronov and Mirlin.^^ 

For a charged fluid it is easiest and most convenient to evaluate Qc in the form 
(10.20): 



ge = -r(i + A,-ir)-i = -(i + rA-i)-ir , 



A ~ A + P 



(11.12) 
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Here is the propagator for electromagnetic fields. In the two-dimensional ap- 
proximation it is given by (10.17). Its form in the ultra-thin film approximation 
can be deduced from (5.1). 



= < 



-3D 




q'^ — 



(ultra-thin film approx.) 



(two-dimensional approx.) 



(11.13) 







Here esD is the coupling constant in three dimensions. 

We evaluate Qc in the two-dimensional approximation. We introduce 



B = 



27r 1 
'Nq 



Then 



-A 





(11.14) 



(11.15) 



The substitution of (11.5) and (11.15) into (11.12) immediately confirms that 
the a-dependence entirely drops in (or Q^), since both {T°'',T'^'') and {T^'^,Tl''^) 
are proportional to (g,w). It is a consequence of the current conservation (11.2). 

The computation of Qc is rather involved. We record intermediate steps for 
readers' convenience. Ar^F is found to be 



I q{qAno - BHi) u){qAno - BHi) i{qAni - SHa) \ 

-wgAHo -ijo'^AHq -iLoAHi 

\ iq{qBUo + AUi) iLo{qBUo + AUi) -qBUi - Alia / 
Further the determinant of (1 + Aj^F) is evaluated to be 
Ac = det (l + A-^F) 

= l + {q'^ - u'^)Ano - 2gBni - AU2 



(11.16) 



+ {q'{A' + B') - wM^KHf - H0H2) 



(11.17) 



= (^)'(ni -non2) + e2Ho 



— up' 



n2 + 



r(n2-non2) 
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Straightforward manipulations lead to 



(1 + A-ir)-i = — ■{! + Ml + M2) 





Ml = I 

iq'^B -lujqB {q^ - uj'^)A ^ 

-qB -iqA^ 

-qB +iu)A I Hi 

—iqA —itoA —qB 



-A iB 

-A 




(11.18) 



-w^ -ojq 



M2 = 



ujq 




q' I AUo 




a; 
-g I Bill 





LoqA'^ q^{A^ + B^) -iuqAB \ (H? - noHa) 
/ 



In computing F (l+A^ ^F) ^ , the contribution from the M2 part vanishes thanks 
to the current conservation F M2 = 0. The final result takes a simple form. 



Qc = 



1 

a: 



q^A ujqA 


ojqA uj 




iq^B iojqB 












ujq 


q^ — up- 




27r i 











\iq 



ujq 













)-: 





• 










-iq 
-iu! 



{Ill - II0U2) 



(11.19) 



It is easy to check that the formula (11.19) reduces to (11.11) in the neutral 
limit 0. Indeed one sees that 



Ae ^ A„ = ( — j {Ut - U0U2) 

Hi - ^ {Ul - U0U2) = -Hi - ^ A„ 



(11.20) 



As we shall see, the difference between Q„ and Qc is important in discussing super- 
conductivity. 
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Ho 


= r"°(c.,q) 


= r(^)°"(c.,q) 


iglli 


= r"2(a',q) 


= rW"^(-,q) 


Ha 


= r2^(a',q) 





12. Evaluation of the kernel 

We evaluate the kernel T'^''{ui, q) defined in (10.7), or equivalently the invariant 
functions Ila's in (11.4). Wc need to evaluate the four diagrams (b), (c), (d), and 
(e) in Fig. 8. For q = (g, 0) we have 



(12.1) 



The diagram (b) has been already evaluated in (10.12): 

T^^^^''(u;,q) = -—6^'' . (12.2) 
m 

To evaluate other diagrams, we first examine the zeroth order propagator 

G{x,y) = -i{T[^{x)^^y)]) 

a 1 -2 \ , (12.3) 

Here the expectation value is taken with the ground state (6.13) and ipi^) satisfies 
the mean field equation. 
More explicitly 

Gix,y) = -ifeixo - yo) ^ -0{yo -xo)J2 ] «a(x)«a(y)*e-'^"(^«-'"') (12.4) 
where Ua(x) is given by either (6.6) or (6.8). In the Landau gauge 

|JV|-1 



G{x, y) = -il e{xo - yo) ^ -9{yo - xo) ^ > 

n=\N\ n=0 ^ 



n=\N\ n=0 

^^,-2.ip(xi-j;0/ii / ^2-a^2 ^„. f y2-y2 \ (12.5) 
P 

= e"^(^'«) • Go{x - y) . 

Here 

(j){x, y) = -e{N) ^ {xi - yi){x2 + 2/2) 



f 00 \N\-1 . 

Go{x) = -il e{xo) Yl -^(-^0) Yl \ 

^ ^ — I Ar\ n — n 



X 



1 

" J -00 



(12.6) 



-z{x2) = e{N) II 



41 



Yutaka Hosotani 



We have introduced a new integration vbariable z = {2-Kpl/Li) — e{N){x2 + y2)/2l, 
taking the limit Li oo. It should be noticed that the Green's function G{x,y) is 
not manifestly translation invariant, but is invariant up to a gauge transformation, 
due to the presence of a non-vanishing magnetic field. The Fourier transform of 
Go{x) = Go{t,x) is given by 

Go(p) = Go(a;,p) = J dtd-K Go(i, x) e*('^*-P'') 



where id{n) = 



+ie for n > \N\ 
-ie for n < \N\ 



We come back to evaluating the remaining diagrams. Recalling in (10.4), 
we see that the diagram (c) in Fig. 8 yields 



2i 



dxdy G{x,y) [-aoiy)]G(y,x) [-ao(x)] 
= ^ j dxdy ao{x) Go{x - y) Go{y - x) ao{y) 



2i 



- I dxdy [-ao{x)] -—a^^^^y) 



2m 



X ( Gix, y) ■ [5| - ia^ (y)]G(y, x) - [aj + id^ {y)] G{x, y) ■ G{y, x) ) 

With the aid of 

[aj - id^x)] G{x, y) = +e'^(-'y) DjGo{x - y) 
[aj + id^iy)] Gix, y) = -e">'^-'y) DjG^{x - y) 

DfGoix -y)= (9J T ie{N) e^" ^^^) Go{x - y) , 
eq. (12.10) becomes 

~ 4m / ^^^^ cio{x) a^^'>' (y) 

X [Go{x - y) ■ DjGoiv - x) + DfGo{x - y) ■ Go{y - x)} . 
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(12.8) 



so that 

rWoO(g) =ij^ Go(p) Goip - q) . (12.9) 
The diagram (d) in Fig. 8 yields 
i j dxdyaoix)T^''^''^{x,y)a^'^^{y) 



(12.10) 



(12.11) 



(12.12) 
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Note that T'^'^^°^x,y) is a function of x — y only. Hence, in the Fourier space 

^^'^°'^^) = ~2^/(S^ {Go{p)-DTGo{p-q) + D+Go{p)-Go{p-q)} (12.13) 
Similar manipulations lead to 
i J dxdy a(i)J (a;)r(^)^'^(x,y) a^^^''{y) 

{i^fc Go(y - x) • DjGoix - y) + D+Go{y - x) ■ D+Go{x - y) 
+DfD^Goiy - x) ■ Go{x - y) + Go{y - x) ■ DjD+Go{x - y)} , 



(12.14) 



and therefore 



i f d^p 

"i^y (2^ 



{d-Go{p) ■ DjGoip - 9) + D+Go{p) ■ D+Go{p - q) 
+D+D^Go{p) ■ Go{p -q) + Goip) ■ DTD+Goip - q)} • 



(12.15) 



We need to evaluate r^''^°°{q), T'^'^^°^{q), and T'^"^'^'^ {q) in the frame q = {q,0) 

to find Ho, Hi, and 112. Upon inserting (12.7) into (12.9), (12.13), or (12.15), we 

encounter an w'-integral 

— r duj' 



27r 



[uj' - e„ + iS{n)] [uj' - lo - €m + iS{m) 



= < 



1 



Cn — Cm — OJ — te 
1 

Em - En + W - ie 



for n > |7V|,m < \N\ 
for n < \N\,m> \N\ 

otherwise. 



(12.16) 



Further, since 



DtGoip) = {ip2 ± e{N) A) Go{p) 



n=0 



w — e„ + iS{n) J I 



dX2 



iN)^-^){vnipil + z)Vn{pil-z)} , 



43 



Yutaka Hosotani 



one finds that 



( (1) (12-17) 

^ I +Vk ' {pil + Z) Vn{pil - Z) 
I -Vnipil + z) Vn\pil - z) 

where Vn^ {z) is the p-derivative of Vn{z). Another integral which frequently appears 
is 

C^i{a) = r dx v^^\x) v^{x -a) = {-If H dx v^{x) viP\x - a) . (12.18) 

J — oo J — oo 

It satisfies that 

We start with computing IIq: 

,^n„ = rM»°{.^,5,0) 

= */^ £ „,-,„Va(„) /f e— ..b.i + .-WI«nM-^Wl 

n=0 



CO 

X 

TO=0 



^ . OO \N\-1 oo \N\-1 

+ E E 



n=|JV| m=0 m=|JV| n=0 



(12.20) 

Therefore 

, oo |JV|-1 

^^no = ^ E E 

„=\N\ m=0 (12.21) 

I — ^ + ^ — ^1 c^iiQir ■ 

Evaluation of Hi proceeds similarly. In view of of (12.13) and (12.17), one needs 
to make small modifications to (12.21) 

X 

2m I (^2.22) 
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to find 

oo \N\-1 

2-KmP 



n=\N\ m=0 (12.23) 

To find 112, we first note that partial integrations in (12.15) lead to 

V^^^^\oj,q,Q)=-^ j ^ D^G,{p)-D^G,{p-q) . (12.24) 
This time we modify (12.21) such that 

V mVV Z y ^^2.25) 
Since Ha = r(^)22(^^ q) + r(^)22(^^ g), we find that 

U2 = 

m 

^ oo \N\-1 

n=|JV| m=0 

{ + ^7 -} C^lillf ■ 

(12.21), (12.23), and (12.26) are the results in the RPA and linearized SCF. 
These H^'s, through the formulas (10.15), (10.21), (11.11), and (11.19), determine 
the response to external perturbations. As is evident from the discussions in Sections 
10 and 11, it describes a response to harmonic perturbations aext{x) oc g'P^"*'^* in- 
troduced from t = — oo tot = +oo. One may introduce a perturbation adiabatically 
from t = — oo to the present (but not to t = +oo). In this case the response function 
is related to the retarded, but not time-ordered. Green's function of currents. This 
amounts to making a change^^'^'' 



e„- em-oj -ie Cn - Cm + u! - ie 



(12.27) 



En - Em - w - ie + w -|- ie 

in all formulas. Hence, for instance, 

|JV|-1 



n=|iV| m=0 (12.28) 
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The a;-integral in Cnm{a), (12.18), can be done to yield 

^™'"^") = i""^""^^l^^""^"^) (12-29) 
(a) = ^ { - C(°^i,„(a) + C(°i,,„(a)} etc. 

Z/"(a;) is the Laguerre polynomial. This expression is useful to investigate the 
response function at finite ql. 

For a small momentum ql <^\, one can expand Vm{x — ql) in a Taylor series to 
find 



dxvn{,x)v^j^\x) , (12.30) 



°° 1 /-C 

fe=o 

where the coefficients djj^'s are given by 
dP =5 

dnm = {VmSn,m-l " V m + 1 6n,m+l} 

dnm = ^|\/TO(m- l)5n,m-2 - {2m + l)6n,m + A/(m + l)(m + 2) 5„,m+2| 

c^^m = I y/m{m- l)(m- 2) ^„,m_3 - 3m^/^ (5„,m_ i 

+ 3(m + 1)3/2 _ + + 2)(m + 3) Sn^m+s} 

dnm = ^|-\/™---(™-3)<5n,m-4 - 2(2to - l)-v/m(m- l)(5„,m-2 

+ 3(2m2 + 2m + 1) - 2(2m + 3)^/{m+l){m + 2)Sn,m+2 

+ VTmTTP"^"(mT4) 5„,™+4| 

rf^m = I v'm • • • (to - 4) 5„,m-5 - 5(to - l)'''/2 v/,„(,„, _ 2) (5„,„_3 



+ 5(2to2 + l)y/mSn.:n-i - 5 (2to2 + 4to + 3) /»7TT f5„ . ,„ _i 

+ 5(to + 2)3/2-v/(™ + l)(™ + 3) V(m + l)---(TO + 5) (5„,™+5 } 

(12.31) 
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In particular, for n > m, 

^nmW 2 "''"+1 \ ^Ori,m+lH 0„^m+2 j - 

Qm(«) Qm(«) = a 2 + " I Y + g ')n,m+2| + 

4fn(37n2 + 5) n(n - l)(2n - 1) 

+ I On,m+l 0„^m+2 



96 12 



, rz(n-l)(n-2) 

"1 On,m+3 J + 

(12.32) 

Applying (12.32) in (12.21), (12.23), and (12.26), one finds that for small fre- 
quency and momentum 



2 

N\ 1 27mp 



^2 



(12.33) 



27rne mZ^ \N\m 



We also note that 



1 , o_ = , N' 



2 



Ha = ^ {w^Uo - Hs) = -— + •• • . (12.34) 

gr^ 2-Km 

There are two notable cancellations in the above formulas. First, in IIi, the 

dominant term is exactly N/2Tr so that IIi = IIi — {N/2n) vanishes at g = a; = 0. 

This fact is phrased in the litcratnirc that the bare Chcrn-Simons term is exactly 

cancelled by the one-loop correction. Secondly, 112 vanishes at g = a; = 0. In 
other words, the first term (diagram (b)) in (12.26) is cancelled by the second term 
(diagram (e)). We have mentioned about it just below eq. (10.12). 

n(a;, g = 0)'s can be evaluated in a closed form. It is straightforward to find 



llo{uj,q = 0) 
ni(w,(7 = 0) 



277 W2_^2 

N 



27r w2_^2 

N .2 ^''-''^ 



2tt ul-uj^ 
\N\ uj^iJ^ 
27r a;?-a;2 • 
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In the literature these kernels appear in different notations. For the sake of 
readers' convenience we have summarized them in Table 1 below. 

Table 1. The comparison of various references. Relations among life's, 
r, Qn, and Qc are given by (11.4), (11.11), and (11.19). At 
T ^ 0, the notation Ilj^ and has been adopted in this 

artick; and in rt;f. 38. 



This 


ref. 27 ref. 29 ref. 49 


ref. 32 ref. 44 


ref. 38 


article 


FHL CWWH FH 


HC HHL 


RSS 


T = T^O 


T = T = T^O 


T = T^O 


T = Tt^O 




uj^O w^O uj^O 


w = a; = 


a; ^ a; = 


no 


_ (^Y^y: 

\2ttJ Ue 




Ho 


Hi 






Hi 




--(I + S2) 
m 






Ha 






-n2 










or 




e-'i?r 




or 









13. Phonons and plasmons 

It follows from (10.11) or (10.21) that the location of poles of the response 

function Qn or Qc in the Fourier space determines the spectrum of excitations 
which couple to currents J^. In SCF they appear as self-consistent configurations 
in the absence of external fields. Indeed, 

Qn ^ Jind = dext = 

or (13.1) 

Qc ^ Jind = C^ext = 

has a non-trivial solution (Jind 0) only at (Wjq) for which detQ~^ = or 
detQ-'^ = 0. 
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From (11.11), (11.17), (11.19), and (11.20) one finds that the location of poles 
are determined by 

An = (^)'(ni - noHs) = (13.2) 

for a neutral anyon fluid, and by 

/27r\2 _9 _ „ (=2 — „ — 

= ( 1^ ) - noHa) + e^Ho - ^ + ^ 2 (^i - ^0^2) = (13.3) 

for a charged anyon fluid. 

Solving (13.2) for a small momentum, Fetter, Hanna, and Laughlin first showed^^ 
that a neutral anyon fluid admits a phonon excitation and is a supcrfluid. Later 
the equation was numerically solved for a flnite momentum in refs. (51) and (55). 
The spectrum in the charged case has been examined in ref. (58). 

In this article we confine ourselves to the spectra at small momenta <C 1. For 
a neutral fiuid one can employ the formula (12.33), as is justified a posteriori. Since 

Ho = 0(1) 

Ti^=li^-^ = 0{io\q^) 
IL2 = 0{uj\q^) , 

Eq. (13.2) is solved for a small momentum by 1X2 = 0, or 



2 22 
u =c,q 



Cs = V iv ujcl = n . 

m 

We have recovered Ti in the last relation. It is a phonon excitation. The velocity 
does not depend on A'' in this approximation (RPA, linearized SCF). 

In the charged case u approaches a finite value as ^ 0. To find Lo{q = 0), we 
insert (12.35) into (13.3). Many cancellations take place. One finds, at g = 0, 

2^^^n?-non2)= 



N 



e'Uo + —Ii2 =2 22 
2-K — uj^ 

-^(n?-non2) ' 



^ ^ V 27r / a;2(a;2-a;2) 

Hence Eq. (13.3) reduces to a polynomial equation for a;^: 



2 I Arl„2 \ 2 

= . (13.5) 



o/_ _ \N\e 
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In other words, the dispersion relation is 

[\N\^ le^Ue 
V ZTT V m 

This is nothing but a plasmon, representing a plasma oscillation. There is only one 
solution for u. [The two solutions in ref. (38) are the result of 11^(0;, 0)'s being 
approximated by 11^(0, 0)'s.] 



14. Hydrodynamic description 

In the preceeding sections we have integrated the matter field first, to obtaine the 
effective theory for the gauge fields. The kernel F'^'' which appears in the effective 
action (10.2) or (10.6) has an important physical meaning. It appears in (10.14) as 
the coefficient relating gauge field configurations to the induced current. In SCF, 
for a given gauge field configuration a^(a;) (in the neutral case), 

j^x) = irix)) 

n f n^ (14.1) 

= neS>'°- J dyT^''{x-y)a';}\y) + --- , 

with which the field equation 

-^e'^'''f.p = J''{x) (14.2) 

has been solved sclf-consistently. 

One can read cq. (14.2) differently. It says that the Chern-Simons field strengths 
are nothing but the currents. The roles of the two equations (14.1) and (14.2) are 
interchanged. Now the latter gives a relation, while the former yields equations for 
the currents J'^{x). 

The time component is the density field J^{x) = n{x), while the spatial com- 
ponents define the velocity fields v(a;) by J*^ = n?;*. The equations for n{x) and 
v(a;) give the hydrodynamic description of the system. This is the viewpoint origi- 
nally adopted by Wen and Zee.^" The advantage of this approach is that everything 
is expressed in terms of macroscopic physical quantities so that physics is clear- 
est. In particular, as wc shall sec, it gives an interpretation of a phonon excitation 
as a breathing mode of a density wave, the picture first spelled out by Wen and 
Zee. Anyon fluids are unique systems in which the "microscopic" RPA or SCF is 
equivalent to the "macroscopic" hydrodynamic description. 

Substitution of (11.4) into eq. (14.1) yields, in the linear approximation which 
drops 0[(Jind)^], 

J° =ne-%/ofcno + 6(i)ni , 

(14.3) 

7*= = -iLO fok no + e'^^foi Hi - ie'^^qi b^'^ , 
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where b^^^ = —f[l^- At T = all n(w, g)'s have finite limits at w = g' = so that the 
right sides of (14.3) are expressed solely in terms of the field strengths Hence, 
with the aid of the fundamental identity (14.2), cq. (14.1) gives a set of differential 
equations for J'^{x). We remark that at finite temperature Uq develops a pole 
so that eq. (14.1) becomes an integral-differential equation. We shall come back to 
this point later. 

Eq. (14.3) becomes 

-Hi e^^Jt^ = -Ho iuj 4^ + Il2 iqj JP ^ 

where Hi = Hi - (iV/27r). A solution to (14.4) exsits only if Hi (fii - IL0II2) = 
where 112 = w^Hq — q^Tl2. Hi = is not permissible, since the equations would 
imply 112 = as well, which is incompatible. Hence we have 

Hi - HoHa = . (14.5) 

This is the same as eq. (13.2), A„ = 0. As was shown in the previous section, 
it admits a phonon spectrum (13.4). We also note that eq. (14.4) contains the 
continuity equation. Indeed it follows from (14.4) that 

Hi (wJ^d - ft^d) = . 

Here the fact Hi 7^ (at finite oj and q) is important. 

To recognize that the phonon mode found in (13.4) represents a density wave, 
we write the equations in terms of n{x) and v(a;). In the long wavelength limit, 

one may insert (12.33) into (14.4). Since Hq, Ha =0(1), while Hi = 0(0;^), the left 
sides of eq. (14.4) give negligible contributions. Hence we have 

doJU^)+cld^jU^) = 

with the continuity equation 

doJ?^di^)+dkJUx) = . (14.7) 

The first of (14.6) says that there is no circulation (vorticity). The second of (14.6) 
and (14.7) give 

{d',-cldl)jPM = o , 

implying a density wave. 

It is easy to see that (14.6) results by linearizing the hydrodynamic equation. 
The Euler equation for an ideal fluid (with no viscosity and thermal conductivity) 
is 

dv 1 

— + (v-V)v = VP (14.8) 

ot mn 
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where P{x) is the pressure 



KdVJr \dV)s 



(14.9) 



We have evaluated the energy density E/V at zero temperature in Sections 6 
and 7. In the RPA and SCF 



SCF/RPA : 



E = V 



TTiVf 



m mV 

where A^e is the total anyon number, and in the Hartree-Fock approximation 

1 ttA/-? 



(14.10) 



HF : 



E 



2 mV 
29 ttN^ 
32~mV 



Hence 



or 



7rn 
m 



1 TTn'^ 

2 TO 

29 Trn^ 



for N = ±1 



for N = ±1 



for N = ±2 



in RPA/SCF 



in Hartree-Fock 



(14.11) 



(14.12) 



(14.13) 



for A^ = ±2 

32 TO, 

Substituting (14.12) into (14.8) and keeping only linear terms in n and v, one 
finds 

dv 2tt „ 
■ Vn 



dt 



do J" = doinv'') 



2Tme 



(14.14) 



which is exactly the second equation in (14.6) with Cg given by (13.4). This deriva- 
tion demonstrates that in the Hartree-Fock approximation the sound velocity is 
modified to^^ 

for A = ±1 



m 



(14.15) 



29 flJ-KUe 



for A = ±2 



4 TO 

Crucial in the above argument is the fact that the energy density is given by 
(14.10) or (14.11), independent of the number density n. The energy density is 
proportional to n^. In the Hartrcc-Fock language |A| lowest Landau levels are 
completely filled even for slowly varying density n{x). As n{x) varies, the magnetic 
length I also varies such that the particles precisely fill the space. If one looks at 
the motion of each particle (in the half classical picture) , the Larmor orbit expands 
and shrinks periodically as the density changes. It breathes. Wen and Zee called it 
the breathing mode.^'^ 
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15. Effective theory 

The effective theory of anyon fluids in terms of Chern-Simons and Maxwell fields 
is obtained by integrating the matter field ^l>{x) first. We have already encountered 
it in Sections 9 and 10. It is given by 



Aff (a, Af=" = £^^{a) + £^^\A) + £f (a + eA) 



-0 -^0 

where C^^{a) and C^^{A) are defined in (9.1): 



(15.1) 



C^^[A] 



477 



4 A"^ 



eneAo . 



Lf{ci + eA) summarizes the effect of the matter field, and is given by (10.6) with a 
replaced by a + 6^4 : 

CF{a + eA) =6p {a + eA) + cf {a + eA) + --- 



Cp^ (a + eA) — —He (ao + eAo) 

{a + eA) = +\ (flW + eA^) T''^ (a^i) + eA,) . 



(15.2) 



In the linear approximation of SCF, or equivalently in RPA, higher order terms in 
a^^^ + eA are neglected. 

The kernel T^" has been evaluated in Section 12. If one is interested in physics 
at a large length scale, F'"' can be expanded in a Taylor series in 9^. In this section 
we shall retain only the most dominant terms. From (12.33) and (12.34) 



no = 

Hi = 

n2 = 



47r^ne 
N 
2^ 

27rm 



for 



\n\qI = 



\N\ 



\/27rne 

|iV|m 
27rne 



g < 1 



•w < 1 



(15.3) 



We note that the linear approximation gets better for a larger |iV|, but the low 
energy approximation (15.3) breaks down when \ 'N\ becomes too large. 
Since 

Cc^{a) = rie ao - ^ e^^P a^^:^ d,Ap , 

all linear terms in C,f.-s{a,-A) cancell each other. (11.4) and (15.3) immediately lead 
to 

£eff (a, A)^=« = -l-e^^f>afd.af - i F^, + (a + eA) 



6l\a + eA) 



47r 

-— e^^>\a^^^ + eA^) + e^p) 



Am 



(15.4) 



+ (/o,- + ei^o.)^ - (6W + ei?)^ 
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Or, noticing the cancellation of the bare Chern-Simons term, one may write 

(15.4) or (15.5) is the effective theory of a charged anyon fluid, valid for slowly 
varying configurations. It replaces the Ginzburg-Landau (GL) free energy for BCS 
superconductors. Instead of the GL order parameter \E'GL(a;) we have Chern-Simons 
gauge fields a^^^ {x). Higher order terms, namely terms cubic or quartic in a^^^ -|-eA, 
become important for large gauge field configurations such as vortices, but have not 
been evaluated so far. 

The effective theory (15.4) or (15.5) was first derived by Hosotani and 
Chakravarty^^ for static configurations at T = 0. There is an alternative way 
of writing an effective theory. Introducing a scalar field <j){x) in place of the Chern- 
Simons field a^(a;), Chen et al. have written down^^ 

£eff(0, A)«™ ^ _ 1 ^2^ ^ g^p^^^^ _ 

+ \{^ - 92Aof -\cl {djcl> - g2Ajf (15.6) 

Ne I m 

with the soimd velocity Cg defined in (13.4). It has been known that in spite of 
different forms both (15.5) and (15.6) lead to the same predictions for many physical 
quantities. Similar effective Lagrangians have been written by Pradkin^^ and by 
Banks and Lykken.^^ 

As we shall see in later sections, the effective theory £eff (a, A)^"° can be easily 
generalized to finite temperature. There we shall find that not only the coefficients 

(2) 

in Cp {a + eA) become T-dependent, but also a new important term proportional 
to (flo + eAo)'^ will turn up. 

We notice that the effective theory £gff{a,A)'^^'^ neatly summarizes the self- 
consistent field (SCF) method. Equations derived by taking variations over a|/^ (a;) 
and Ai^{x) are 

^H. ft^i^p f(i) — r/* 

4^ £ Jiyp - -^ind 

(15.7) 

^L(^) = -7W- {4'^(« + e^) + ---} 

6a}i'{x) 

= -r'*'' (a(,i) + eA,) + 0[(a(i) + eA)^] 
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The last equation which expresses the induced current J/^^ in terms of the two 
gauge fields and may be viewed as a new London equation. The first and 
second equations lead to an identity 



_g^^.p_^(i)^3^^.M, (15- 



or in the component form 



In ^1^-9) 

— fok = e'''doEi+dkB , 

with which Chern-Simons fields may be eliminated. We note that the identity (15.8) 
or (15.9) is valid beyond the linear approximation. It is indeed a direct consequence 
resulting from the general structure of (15.1), the minimal gauge couplings. 
For slowly varying configurations 



JlM = +^ [b^'^ + eB) - ^ a, (/o, + eF,,) + • • • 

N N-^ 
JUx) = +^ e"' {foi + eFoi) - ^ e'^'diib^'^ + eB) + ■ ■ 

With the aid of (15.9) we eliminate Chern-Simons fields to obtain 



(15.10) 



(15.11) 



eJ-nd = { - doEk + e'^'diB } - |— e'^'OiB 



27rm 

+ ^e'''Ei--€'''di{div-E) + --- 

ZTT m 

where rotE = diE2~'d2Ei. Notice that the dominant terms in Eq. (15.10) represent 
an integer quantum Hall effect in the system.^^ 
Inserting (15.11) into the Maxwell equations 

divF, = eJ:°H 
-doEk + e''' diB = eJ.^^ , 



one recognizes that the terms in parenthesis { } in (15.11) exactly cancell the left 
sides of the Maxwell equations. Crucial in this cancellation is the fact that the 
coefiicient of the induced Chern-Simons term for a^^^ + eA^ in (15.4) is exactly the 
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negative of the coefficient of the bare Chem-Simons term for a^. We thus arrive at 
equations 

TO „,\ „ Nm ,. „ TO „ , 
1 - ^— B - div E - ^— do (roi E = 



2-iT N 

Ek-^dk (divE) dkB = 

e^TO TO 



(15.13) 



which describe electromagnetic fields in anyon fluids. 

16. Meissner effect at T = 

We examine a response of a charged anyon fluid to an external magnetic field. 
If the anyon fluid is a superconductor, a sufficiently small magneic fleld must be 
expelled from the system. It must have a Meissner effect. 

We shall show that it is indeed the case at least at T = within RPA and SCF. 
There are two ways to demonstrate it, one in the real conflguration space with an 
external magnetic field applied outside the body, and the other by introducing a 
test current of a ^-function type in the middle of the body. The former approach 
corresponds to solving the Ginzburg-Landau equation in the BCS superconductors, 
whereas the latter to a linear response theory.^^ We discuss both. 

First we suppose that an anyon fiiud occupies a half plane (xi > 0, — oo < X2 < 
+oo) and an external magnetic field is applied such that = i?(x) = Bgxt for 
xi <0. The problem is to determine the magnetic field configuration -B(x) = B{x-i) 
inside the anyon fluid. One expects damping behaviour B{xi) oc exp(— a;i/A) if the 
anyon fluid is a superconductor. 

In this section we consider a "sufficiently small" B^xt- We expect that deviations 
of both Chern-Simons and Maxwell flelds from the ground state values are small 
so that Eq. (15.13) may be employed. Together with the boiuidary condition 
B{0) — i?oxt, the magnetic field B{xi) (xi > 0) inside the system is determined. 

For configurations under consideration, Eq. (15.13) becomes 

e^TO TO 

E2=0 . 

At this point one has to examine numerical values of various parameter. We give 
a summary of numerical values in Sections 18 and 20. It will be seen that to very 
good accuracy (16.1) is approximated by 

E^^O (16.2) 
E2 = . 
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Hence the solution is 

B{xi) = Sext 6-"=^/^ for a;i > 
e^ne 



(16.3) 



The magnetic field is exponentially damped from the surface. The charged anyon 
fluid exhibits a Meissner effect at T = 0. The penetration depth coincides with the 
London penetration depth in BCS superconductors. The persistent current flows 
along the boundary. 

Next wc shall examine the same problem in the linear response theory. We 
imagine that a charged anyon fluid occupies the entire space in the X1-X2 plane. We 
introduce an external current of a (5-function type at xi = 0: 



(16.4) 



which generate an external magnetic field 

Be^t{x)=Boe{xi) . (16.5) 

In the momentum space 

ej|,t(c.s q) = -2Bo ■ (2^)2 5(a;)J(g2) 

-Bext(w, q) = ^ ej2^t(w, q) ^^^ ^-^ 

Al^ti^, q) = 4 e Jext('^. q) > ^e^t = ^ext = • 
1l 

The response of the system to an external perturbation is described by the 
response function Qc determined in the preceeding sections. The relation to the 
induced current is given by (10.21): 

J^/"'^^'^ = -r'*-(aW +eA,) =Q^^ eAr* . (16.7) 
For the configuration (16.6) 

7.° , - - p ^2 _ _ _ n02 72 

1l 

JL = (16.8) 



t2 __r)22p/l2 _ _ f_ /2 
■^ind ~ '^^ext ~ „2 "^ext 

1l 
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At this point we need to evaluate Qc'^ for a; = 0, ^2 = 0, and small qi = q. We 
notice that 

no,ni,n2 = o(i) 
ni,n2 = o(g2) 

n2-non2 = o(i) 
Hi - non2 = 0(g2) . 

Hence in Ac in (11.17), only the second and fourth terms are relevant. Explicitly 

(16.9) 

fN\2 -2 2, 

In the second line we have suppressed numerically negligible terms. A is given in 
(16.3). 

Prom (11.19) 

Qf = ^{-n2 + e^(nf-non2)} 
Qf = ^{-n. + |(n?-non2)} . 

It is straightforward to see that in our approximation 

n'22 -i 



1 

1 + aV 



iN g5 
4e^l + AV 



(16.10) 



The total current in the presence of the perturbation (16.4) becomes 



72 _ t2 , t2 

e2 



1 2 ) "^xt 

1l ' (16.11) 

2 

9? j2 



-2 „ e'^t ■ 
1 + A (?2 

Notice that J^^^ vanishes at q\ = 0, i.e. the external current is completely shielded. 
The total magnetic field is given by 



5tot(w,q) = — eJ^t 
91 



2 (16.12) 



A <7i 



= -2zi3o ^ • (27r)^ (5(u;)(5((Z2) 



1 + A 



2 
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SO that in the configuration space 

=Boe(a;i)e-l^^l/^ . (16.13) 

We have reproduced the same result as in (16.3). The Meissner effect is complete 
at T = for a sufficiently small external field. 

For completeness we look at a response to an external static charge: J^^t 7^ 0, 
Jgxt = 0. One finds that, for w = 0, 

Qc° = -7^{9'no + e2(n2-non2)}~-4 • (i6.i4) 

Noticing 
we find that 

J°„, = eA°,, = ^ g°° J°,, ~ - J°,, . (16.15) 
An external charge is completely shielded, as it should. 

17. T 7^ — homogeneous fields 

The behaviour of anyon fluids at finite temperature is particularly interesting to 
know. In a sense the behaviour of anyon fluids at zero temperature is very similar 
to that of conventional superfluids or superconductors. For instance, a charged 
anyon fluid exhibits a Meissner effect for sufficiently small magnetic fields with the 
same penetration depth as in BCS superconductors. It is very difficult to sec an 
effect of the unique structure of the ground state, namely the complete filling of 
Landau levels in the Hartree-Fock approximation. (There is a tiny effect of P- and 
T-violation, which we shall briefly touch on in Section 23.) 

What would happen at finite temperature? Does a charged anyon fluid behave 
quite differently from BCS superconductors? Is a charged anyon fluid really a "high" 

superconductor? We shall show in this and following sections that there is an 
indication that a charged anyon fluid has Tc around 100 K, but not around 5 K or 
1000 K. After all the most important feature of observed high superconductors 
(cuprate superconductors) is that they have Tc ~ 50 — 100 K. 

How is the order of a charged anyon fluid destroyed? In this respect it is suit- 
able to consider effects of both finite magnetic fields and temperature.^'* Recall 
that at zero temperature with no external fields Landau levels with respect to the 
Chern-Simons magnetic field are completely filled in the Hartree-Fock approxima- 
tion. Particles, or holons, feel only the sum of Chern-Simons and Maxwell gauge 
fields. They interact with the gauge fields in the combination of + eA^. 

If an external uniform magnetic field i?ext is applied in the direction perpendic- 
ular to the two-dimensional plane, to the first approximation, particles feel the total 
magnetic field 6tot = b^^^ + eSgxti with which the Landau levels are not completely 
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filled any more. If the relative sign between b^^^ and eSext is negative, less states are 
available per Landau level so that some perticles must be put in the higher energy 
level. On the other hand, if the relative sign is positive, more states are available so 
that there appear vacant states in the top filled level. Hence, if the complete filling 
were essential for superconductivity, an external magnetic field would destroy it. 

Similarly, at finite temperature, the levels are not completely filled due to ther- 
mal excitations. Superconductivity should be destroyed at sufficiently high temper- 
ature. 

The anaysis at finite fields and temperature, however, is complicated by the 
plausible breakdown of the approximation in which homogeneous configurations are 
supposed. It is likely that a finite uniform external magnetic field creates vortices 
in anyon fluids, giving rise to inhomogeneous, though patterned, configurations. At 
finite temperature vortex-antivortex pair creation would become important. 

Effects of vortex formation has been examined by Kitazawa and Murayama^® in 
the case of neutral anyon fluids. They have argued that vortices bring a stability 
to the superfluidity of neutral anyon fluids. At the moment the existence of real 
(electromagnetic) vortices in charged anyon fluids is yet to be established. 

In this and following sections we shall examine effects of finite fields and tem- 
perature, ignoring contributions of vortices. This is a drastic approximation, which 
has to be improved in future. However, we shall find that even in this approxima- 
tion anyon fiuids exhibit interesting behaviour which is quite different from that of 
conventional (type I) superconductors. Some of the behaviour will be modified by 
the incorporation of vortices, but all this is certainly essential for full understanding 
of anyon fiuids. 

The first evaluation of finite temperature effects in the model under consideration 
was given by Randjbar-Dacmi, Salam, and Strathdee.^* Hetrick, Hosotani, and 
Lee^** subsequently confirmed their result, discussing more physical implications 
with additional effects of finite fields. Later many authors, particularly. Fetter and 
Hanna,^® recovered the same result by different methods. 

The evaluation consists of two parts. First thermodynamic quantities are evalu- 
ated for uniform field configurations at finite temperature, from which self-consistent 
uniform fields are determined. Secondly, inhomogeneous deviations of the fields from 
the self-consistent uniform configuration are incorporated in perturbation theory at 
finite temperature. 

With a uniform magnetic fields 6^°^ and B^^\ Landau levels are formed with 
respect to b^^^ + eB^^\ In this section we suppress the superscript to simplify 
the notation. It will be recovered when inhomogeneous configurations are examined 
in Section 19. The number of states per area per Landau level is 

\b + eB\ 



which defines the magnetic length 
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Energy levels are given by 

One of Chern-Simons equations implies that 

, 27rne 

is still valid. Therefore the filling factor as a whole is 



{n = 0,1,2,- ■■)■ 



(17.3) 



(17.4) 



V = — = 27rne • / 
riL 



Nb 



b + eB 



\N\ 



(17.5) 



For \eB/b\ «: 1, 



u=\N\- I 



NeB 
27rne 



(17.6) 



We keep the above definition of z/ at T ^ 0. 

The distribution function pnp at level n with the second index p (defined in 
Section 6) is given by 

where /i and T arc the chemical potential and temperature, respectively. The 
chemical potential is fixed by the condition 



(17.8) 



The summation over p gives V/2'kI'^, as is most easily seen in the Landau gauge. 
Hence 



V = 2itI^ ■ Ue = pn ■ 

n=0 

The energy density £ and entropy density <S are given by 

n,p 

2TTnl / 1\ 
v^m ^ — ' \ 1/ 

ri 

•5 = -:^^ {p„lnp„ + (l-p„)ln(l-p„)} 



(17.9) 



(17.10) 



n,p 



= --^^ {Pnlnp„ + (l-/9n)ln(l-p„)} 
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The free energy density J^[B] is 

J^[B]=£-TS . (17.11) 

A good approximation to p„ is obtained by examining numerical values of various 
parameters. As we shall see in the next section, with typical values of m ~ 2me 
and rie ~ 2 X lO^cm-^, 

e=^ . ^•2800 K 

\N\m \N\ 

1. 27rne 2 ^ ^^'■^^> 

-b=^T7^ ~ ■nr77-1200 T . 

e \N\ \N\ 

With this choice, Tc will turn out about 100 K. In other words, with TV = ±2 and 
for T < 200 K and B < 30 T, the filling factor v is very close to jA'^j and thermal 
excitations are appreciable only to the |A''|-th Landau level: 

f 1 forn < lATI - 2 
~ < (17.13) 
\ forn > |7V| + 1 . 

Hence Eq. (17.9) becomes 

X = u - \N\ = {p\N\_i - 1) + p\N\ . (17.14) 
It is called the two-level approximation. The condition for its validity is 

e""/"^ «: 1 . (17.15) 

Since 

_ 1 _ 1 

Eq. (17.14) is solved by 

^(x,T)±i^e±(l^l^-^)/^ 

T±^ ^T^cosh- 

which determines the chemical potential /x with given T and B. (Note that e is a 
function of S.) 
In particular, 

/i(r) = lA^le at S = . (17.17) 

There is no T-dependence in the two-level approximation. The correction has been 
evaluated by Randjbar-Daemi et al. to be^^ 

^(S = 0) = |Ar|e-|e-l^l^/^ + --- , (17.18) 
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which is exponentially small. 

We also remark that T ^ limit is singular, as can be seen from the presence 
of sinh(e/2T) or cosh(e/2T) in (17.16). At zero temperature for small \x\ 

( (\N\ - i) e for a; < 
H(T = 0) = 1 , (17.19) 

I (1^1 + 5) e for a; > 

which also follows from the consideration of the Fermi level. Thermodynamic quan- 
tities at T = have a singularity at x = 0. We shall come back to this point in the 
next section. 

Let us define ^ 

P=P\N\ = ^e./2T + i ' (17-20) 

in terms of which the energy and entropy density in the two-level approximation 
are given by 

5 = -^ { plnp +(l-p)ln(l-p) (17-21) 
+ {p -x) ln(p -x) + {l-p + x) ln(l -p + x) | 

There are two parameters, x and T. For experimentally available magnetic fields 
we always have \x\ -C 1. However, z{x,T) in (17.16) depends on T sensitively with 
small, finite x. It is easy to see that 

1) |a;| «: 1 , \x\ «: 1 : 

2) \x\ « 1 , \x\ » 1 : (17-22) 
a; > : ^; = 1 e"^/^"^ = 1 P\n\ = x 

x<0 : z = —xe'^l'^'^ p|Ar|_i = 1 + x p|Ar| = . 



Note that with the numerical values (17.12) the boader line defined by \x\ e^/"^^ = 1 
is given by, for = ±2, 



X = 


-6 


xlO- 


ll 


B = 


10-3 G 


T = 


60 K 


X = 


-6 


xlQ- 


-9 


B = 


10-1 G 


T = 


74 K 


X = 


-8.3 


xlO" 


-7 


B = 


12.5 G 


T = 


100 K 


X = 


-6 


xlO" 


-6 


B = 


10+3 G 


T = 


146 K 


X = 


-6 


X lo- 


-3 


B = 


10+5 G 


T = 


279 K 



(17.23) 
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It is interesting that the crossover takes place around 70 - 150 K for moderate 
magnetic fields. However, it should be borne in mind that the above result is in the 
approximation which ignores contributions of vortices. 

For completeness we evaluate the specific heat (per volume) and pressure in the 
two-level approximation. From (17.21) one finds 



It follows from (17.16) and (17.20) that 



\dT)v 2r2 (2e'=/2r + i)(^-ie-^/2T + i) 

(17.25) 

a;sinh(e/2T) 



[l + a;2 sinh^(e/2T)] 



1/2 } • 



The expression for the pressure at finite x is lengthy. The result for a; = is 
simple. 



KdVJT 



(17.26) 

(^i»i'-^)-(i).''(^-W"r^)- 



Insertion of (17.20) with z = 1 shows that the two terms in the last parenthesis 
cancell each other. Hence 



m V 



2 



^ ^1 + ^pJ at x = 0. (17.27) 

There results only a tiny correction. However, this may be an artifact of the uniform 
field approximation. We shall see in Section 21 a sign of instability in a neutral 
anyon fluid. 

18. de Haas - van Alphen effect in SCF 

Charged anyon fluids have the structure of Landau levels, and therefore should 
exhibit a de Haas - van Alphen effect^^ when external magnetic fields are applied. 
Of course, an implicit assumption is that the system remains uniform in the pres- 
ence of uniform fields, which is probably not true even with a modest external field. 
Observed high Tc superconductors arc of type II, i.e. vortices are formed. Never- 
theless, it is worthwhile to examine how the system respond to external fields in 
the uniform field approximation. We shall see that the Meissner effect at T = 
for sufficiently small fields can be understood as a part of a de Haas - van Alphen 
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effect, and that an important departure from the Meissner effect of the BCS type 
results both at modest external fields and at finite temperature.^'* 

In the previous section we have evaluated in the two-lcvcl approximation the 
energy and entropy densities as functions of temperature T and magnetic field B. 
In terms of the free energy dnesity T = £ — TS, the magnetization is given by 

M{T,B) = -^^^ . (18.1) 

B is the total magnetic field (magnetic induction). The relation to an external field 
(thermodynamic field) H is given by 

B = H + M{T,B) (18.2) 

which defines a relation between H and B. 

Let us see first what happens in the T — > limit in the uniform field approxi- 
mation. The mean field energy density is easily computed to be 

,(o,B) = !I!i/i+(l^l-)(Vl^l^^) 

T^nl \eB\ eB 



1 + 



m \ 2'iTne \2'Kne' J (18.3) 

\N\-l<v< \N\ {NeB > 0) , 
\N\<iy< \N\ + 1 {NeB < 0) . 



for 



Notice the appearance of \eB\ in the expression. The energy density has a cusp at 
B = 0, as was first noticed by Chen et al.^^ The magnetization is found to be 



eB>0 MiO,B) = -^^[l-2\Nm\Tl)^} 
eB<0 M(0,i.) = +|^{l + 2|iV|(|iV|Tl)^} 



(18.4) 



It has a discontinuity at B = 0. The magnitude decreases from the value M(0, 0^) 
as \B\ increases. For magnetic fields available in laboratories, we always have 
\eB/27me\ -C 1. Therefore, to good accuracy 

en (eB>0 
M(0,B)=T^ for <^ . (18.5) 



2m 



eB <0 



This is nothing but a de Haas - van Alphen eflfect. (See Fig. 10.) The only 
difference from the standard one is that even in the absence of magnetic field {B = 
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0), we have a Chern-Simons magnetic field {b ^ 0) such that we are at the integer 
fining u = \N\. The magnetization reaches its maximum \e\ne/2m at discontinuous 
points. 



Fig. 10 B vs M as T varies, (a) At T=0. (b) For moderate B as T varies. 

What does this mean? When an external magnetic field H is applied in the 
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direction perpendicular to the plane, the solution to Eq. (18.2) is 

H < -H'^ B = H + H'^, M = +H'^ , 

-H'^ <H< H'^ B = 0, M = -H , (18.6) 

K<H B = H-H',, M = -H',, 

where the critical field is given by 

ff; = !£^ ^ (18.7) 

So long as \H\ < H'^, there is no magnetic field {B = 0) in bulk. It is a Meissner 
effect. However, ii \H\ exceeds H'^, a part of the external magnetic field penetrates 
inside the anyon fluid. (See Fig. 11.) As a consequence of the de Haas - van Alphen 
effect, there is the maximum for \M\. 



Fig. 11 H vs B as T varies. 

Of course, all of these results have been obtained in the uniform approximation. 
In reality the formation of vortices would invalidate the above picture. Then arises 
a question which one is smaller, Hci or H'^. 

So far all quantities are defined in the effective two-dimensional space. We 
have been supposing that three-dimensional material has a layered structure with 
interplanar spacing d 5 A. Two-dimensional quantities (denoted by ( )d=2) are 
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related to three-dimensional quantities (denoted by ( )d=3) by 

/e^N _ 1 _ 1 1 

Bd=2 = \/rf-Bd=3 (18.8) 
ni=^ = dni=^ . 

Hence the Chern-Simons magnetic field h in (17.4) and the critical field H'^ in (18.7) 
are given by 

i-h) =^bd=. = Ll^nr 
\e Jd=3 ed=3 Ny a 



<i=3 



2m 

The hole density nf""^ and nf""^, and the spacing d are measured directly. Typ- 
ical values are 

nf=2 = (l~5)xl0"cm-2 

nf=3 = (2 ~ 10) X 10^^ cm-=^ (18.10) 
d = 5 A = 5 X 10"^cm . 
We also note the conversion formulas: 

rUe = 5.1 X 10^ eV = 2.6 x 10^° cm"^ 

1 G = 1.779 X 10^ cm-2 = 6.903 x 10"^ eV^ (18.11) 
1 K = 8.617 X 10"^ eV = 4.367 cm'^ . 

The effective mass m is not directly mcasiirablc. It can be fixed from observed 
values for the penetration depth at T = or Tc- [T^ is discussed in Section 22.) 
Recalling (16.3), one finds 



KT = Q) = \=\l^:^^,^x\^n^ . (18.12) 

Therefore 



^d=2"'e V ^d=3"'e 



m = Anani-^ X{Of . (18.13) 
If one substitutes the values (18.10) and 

A(0) = 1400 A m= (1.4 ~ 6.9) rUe . (18.14) 

We shall see in Section 22 that 

Tc~100K m = (l~5)me . (18.15) 

(Indeed, the ratio nf^^/m is related to A(0) or Tc.) 
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Upon making use of (18.13), H'^ can be written as 

1 1 



Hi 



2e 47rA(0)2 

In the Ginzburg-Landau theory of conventional superconductors of type II Hd is 
approximately given by^* 

2t:Tic 1 , a „ . „n 

where ^ is a coherence length. If one uses the GL parametrization to high Tc 
superconductors, one typically finds^^ A/^ ~ 100 so that ln(A/^) ~ 4.6. In other 
words, 

roughly i?^ ~ i?ci • (18.18) 
The values for h and H' in three dimensions are 



1 9 r)'*=2 

i6=^— ^ • 1.2xlO^G 

e A/^ 2 X 1014 cm-2 

4 X 10^1 cm-3 m 



(18.19) 



Or, with the aid of (18.16) one has 

--(^)"-- . 

The Chern-Simons magnetic field is huge (~ 1000 T), but H'^ turns out to be 
in a modedst range (-^ 50 G). The rough equality (18.18) suggests that vortices are 
formed in anyon fluids, and that the uniform magnetic field inside the fluid, B in 
(18.6), may represent the average field over the vortex lattice. 

The huge magnitude of h might be related to Hc2 of high Tc superconductors 
at T = 0, which is known to be much larger than 100 T. Related to the huge 
Chern-Simons magnetic field is the energy spacing in the Landau level. 

1 27rnf=2 



\N\m 

= lh — TT-^^ ^ • 2800 K (0.24 eV) . 

\N\ m 2 X 10" cm-2 ^ ' 

Generalization to finite temperature (T ^ 0) is straightforward. Magnetization 
M{T,B) is determined by (18.1) with (17.21). Then, the B vs H relation is found 
from (18.2). We have given the result in Figs. 10 and 11. 
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As is seen from the figures, magnetization alomost vanishes around 100 K. This 
is presumably related to the phenomenon observed in (17.23), and therefore need 
more elaboration by incorpolating vortices. 

19. T ^ — inhomogeneous fields 

The effective action or free energy at finite temperature can be evaluated in 
perturbation theory. One way is to write down the partition function in the path 
integral representation as we did in Section 9 for T = 0. This amounts to rotating 
a time t to an imaginary time t = it (0 < t < {3 = 1/T) with appropriate bound- 
ary conditions on the fields imposed. In other words, we have Matsubara's finite 
temperature Green's functions in place of time-ordered Green's functions at T = 0. 

Most of the arguments in Sections 9 - 11 remain intact, provided that the 
frequency iv is replaced by Matsubara frequency iuin- [ w„ = 2n7rT or (2n + l)nT 
for bosons or fermions, respectively, where n is an integer.] In particular, the 
decomposition (11.4) of the kernel F in terms of II's. and the relation (11.11) or 
(11.19) between II's and the response function Qn or Qc are valid after the Wick 
rotation. However, one should note that all these are for r-ordered Matsubara's 
Green's functions so that to relate them to, for instance, a response function in real 
time, one need to make necessary transformations. There are also studies in the 
real-time formalism of finite temperature Green's functions. 

Again we integrate the ferniion field ip first to obtain the effective theory for 
Chern-Simons and Maxwell fields. We expand the gauge fields around constant 
magnetic fields. In the Landau gauge 

a'=(x) = -6Wa.2+a«'=(x) 

A'=(x) = -B(0)a;2 + A«'=(x) ^ ' 

where 

= ^ (19.2) 
and B^^^ is a constant background magnetic field. According to (17.1) - (17.5) 

e(Ar)(6W+eBW) = i = ^ . (19.3) 
The fermion part of the Hamiltonian is 

He[a + eA] = J { ^{Dk0{Dki^) + (ao + eAo)VV } (19.4) 
where Dk = dk — iia'^ + eA'^). We decompose it into the free and interaction parts: 
He[a + eA] = Ho + H,^t = Ho + Vi+V2 



(19.5) 



/dx^(a«'=)^VV 
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where 

= ao + eAo 



Dk = dk- dkiie{N)-^ , 



X2_ 



The matter part of the free energy with given static (f-independent) gauge field 
configurations is defined by 

e-^^'["+^^l = TV canonical e-'^^^I^+'^^l . (19.6) 

In the previous section we have evaluated the zeroth order free energy: 

C ^ ^ — Tr canonical 6 ^ ^ . (19.7) 

In these formulas the trace is taken over states with a fixed number of particles, 
i.e. over a canonical distribution. One can consider, instead, a grand canonical 
distribution to define the thermodynamic potential fig- 

„-l3ne[a+eA] _ rrv. -|^{He[a+eA]-^lN) 

~ (19.8) 



where 

N 



= V = (TV) (19-9) 
(g) = Tr Qel^i^e-He-nN) 

We have adopted notation N for the number operator to distinguish it from the 
coefficient of the Chern-Simons term, N. 

The original computation of ref. (38) was performed for a grand canonical 
distribution. In view of the rie-dependence of Hq through 6^°^ a perturbation 
theory for a canonical distribution was employed in ref. (44). So long as macroscopic 
physical quantities are concerned, there arises no difference between the two. Even 
at the diagram level there is not much difference except for a minor change in the 
fermion propagator. 

In this article we adopt a perturbation theory for a grand canonical distribution, 
which is summarized in the book of Abrikosov et al..^^ We outline the argument in 
the operator formalism, supplementing expressions in the path integral formalism. 

We define finite temperature Heisenberg field operators by 

d 

-^V'(-r,x) = A^V(t",x) V(0,x)=V(x) 

+ |-^(t,x)=M*V^(t,x) V^(0,x)=^t(x) (19.10) 

OT 

M = Dl + ao + eAo - n 
2m 
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where tp{x) and ^^^i) are the operators in the Schrodinger representation. 

If the Hamiltonian is r-independent, then the equations can be integrated 

as 

Mt, x) = e^^^^-"^) e-^^^'-"*' d 

, , if ^ife = . (19.11) 

V^(r,x) = e^(«=-''^)7^t(x)e-(^^e-^iV) dr 

At this point we make one technical generahzation. We aUow that gauge fields 
a^^^ and A^^^ may depend on r, provided that they are periodic with a period /3: 

aW(/3,x)+eAW(/3,x) =aW(0,x)+eAW(0,x) . (19.12) 

(We assume that the zeroth order parts are r- independent.) This is a technical 
device which enables us to probe dynamical properties, namely time-dependent 
phenomena, of the system at finite temperature, through appropriate analytic con- 
tinuation. 

The transformation matrix for a general r-dependent is defined by 

U{t2,ti) = T, cxp I - ^ ' dr {Heir) - fiN)^ , (ra > n) (19.13) 

where T,- indicates the r-ordering. The specification of the ordering is necessary, 
since [He{T), He{u)] ^ for r-dependent gauge fields. By definition 

f/(r3,r2)f/(r2,ri) = ?7(T3,Ti) (rs > > n) . (19.14) 

Further C/(r) = J7(t, 0) satisfies 

d 



(19.15) 



(19.16) 



^^U{t) = -{H,{t-)-ixN) U{t) 
l.U{T)-' = +U{T)-HH,{T)-t,N) . 

In terms of U{t), Eq. (19.10) is integrated to yield 

V'(r,x) = C/(r)-^ V(x) U{t) 
(r, x) = C/(r) - ^ (x) U{t) . 

Indeed, 

( - ^) U{t)-^ V(x) U{t) = U{t)-' [^(x), H,{t) - p^N] U{t) 

= U{t)-^ M'il){yi) U{t) 
= M tp{T, x) . 

In the last equality we have made use of the fact that the differential operator A4 
commutes with U{t). 
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The next step is to introduce the interaction representation. With the free 
Hamiltonian Hq in (19.5), field operators in the interaction representation are de- 
fined by 

V;i„t(T, x) = e^^^"-"^) V(x) e-^^^"-"*) 
The transformation matrix in the interaction representation is given by 



(19.17) 



S{t) = T, exp I - ^ dr' A„t(T')| 



(19.18) 



= Hint [■0int (r, x) , V int (t, x) ; (a + e A) (r)] 
The fundamental operator identity is 

U{T)=e-^^"°-^^^S{T) . (19.19) 
The diagram method is developed on the basis of (19.19). Let us denote 

g-/3no = xr £-'3(^0-'"*) 

(19.20) 

We define the Euclidean effective action Ie [a + eA] by 

Taking a trace of (19.19), one finds 

lE[a + eA]=l3no-\n{SiP))o . 
The Bloch-De Dominies theorem is applied to (5(/?))q, leading to 

lE[a + eA]=pflo-{{S{(3)),-l} (19.22) 

where the subscript c indicates that only connected diagrams be taken into account. 
For static gauge field configurations, one has 

= ^ Ie = ^0 - ^[{S{P))c - l| for static a + eA . (19.23) 

The path integral representation is obtained for Tr f/(/3) in (19.21) by the stan- 
dard technique: 



[a+eA] ^ J exp I - dr j dx jCe[a + eA] | 

_ _ . 1 „ _ (19.24) 

Ce = ipip + —DlijjDkip + (ao + eAo - At)V'V' 

B.C.: V(/3,x)=-V(0,x) , ^(/3,x) = -V^(0,x) 
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As it stands, this expression is formally obtained from (10.1) by Wick-rotating the 
time axis through 90 degrees and imposing the anti-periodic boundary condition 
on the fermion fields. One should remember however that the rigorous derivation 
follows from (19.21), and that the expression (19.24) is for grand canonical distri- 
butions. 

The free propagator for the fermion field tp is defined by 

G(x,y;Ti - T2) = -(T^[^i'int(x,ri) V^int(y,r2)])o (0 < n,T2 < P) . (19.25) 
It is easy to see 

G(x,y;T) =-G(x,y;r + /3) for - /? < r < . (19.26) 
In the Landau gauge, analogously to (12.5), we find 



X e 



(/?) - 1 for < T < /? 
\pnif}) for-/3<r<0 

= e'^(^'«) • Go(x-y,T) . 



Here 



4){x, y) = -e{N) ^ (xi - yi){x2 + 1/2) 



1 f°° 

C?0(X, r)=Y,^] dz e---^/' Vn[z - ziX2)] Vn[z + Z{X2)] 

-p,)r /pn(/3)-l for 0<T</3 



X e 



~z{x2)=e{N) I . 



(19.28) 



p„(/3) for -/? < r < 



As in the T = case, the Green's function G(x, y; r) is not manifestly translation 
invariant, but is invariant up to a gauge transformation. 
In the Fourier space 

Go(x,r) = ^E / Go(...,p)e--^+-- 

Go(w^,p) = j' dr j dx Go(x, r) e^^'^-^-P'^) 

° (19.29) 

V . I e-'P--- Vn\pil + ^{X2)]Vn\pil-Kx2)] 



(2r + l)7r 
where = 
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We also note that 

(V^(x,T)V(x,r))o = G(x,x;0-) 



1 f°° 1 
= E ^ /_ ''nizf Pn{0) = ^Y.Pn (19-30) 

n °° n 

= rie . 

Similarly 

• Dki' - dI^ • V) = . (19.31) 

(19.30) and (19.31) appear in the first order perturbation. 

As in (10.6), one can expand the Euclidean effective action Ie in a power series 
of a(i) +eA(i): 

lE[a + eA] 

= /3 r^o + j dx Ueiflo + eAo){x) 

- j dxdy\{a + eAf^\xm^{x-y){a + eA)^^){y) + ... 

= /3f2o + j dx ne{ao + eAQ){x) 

\{a + eA)^^{-p)T%\p){a + eAi^\p) + .. 



(19.32) 



where 



= (ti,x) , y=(r2,y) , j dx = j dr J dx 



2-Kr 

A relationship between V'^ and the response function at finite temperature can 
be found easily. Let us denote 

iV-x)= V^(T,x)V'(r,x) 

if- - 1 (19-33) 

3 {t,x) = - — |V(t,x) •Dfe'(/'(r,x) -D^V(t,x) • V(T,x)j . 

Then from (19.13) and (19.16) it follows that 

^^^=-Tr[/(/3,r)/(x)t/(r,0) 
dao{T, x) 

= -Tr U{(3,0)f{r,x) 
= -{f{r,x)).TrU{(3) 
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SO that 



dao(r,x) 
oa'=(r, x) 



(19.34) 



Obviously 



(i'*(r, x)) = (i'*(x)) for static a + eA. (19.35) 
We denote the thermal average of the induced current by 

JC,air,^) = irir,^))-S''°rie . (19.36) 
Making use of (19.32) and (19.34), one finds 

■^';:d (^) = - / ^e{^ - 2/) (« + eA)i'\y) + • • • , 
or in the Fourier space 

Ji:,AoJr,p) = -r^Ei^r,p)ia + eA)^J\ur,p) + --- . (19.37) 

As in the zero temperature case, the self-consistent field approximation (SCF) is 
defined by Eq. (19.37) and the field equations for the gauge fields with the source 
replaced by the thermal average (j''(t, x)). 

Of course, for general r-dependent field configurations, appropriate analytic 
continuation of the field equations is necessary. In this article we mostly restrict 
ourselves to physics for static configurations {tOr = 0), for which field equations are 
the same as those at T = 0. In particular, the response functions Q„ and Qc are 
defined in the same way as at T = 0, by introducing static external fields. In the 
linear approximation, in which higher order terms in (19.37) are neglected, one has 

J;,.^ linear ^g^^^ext Q^^^a^ ■ (19.38) 

The relation among F , Qn, and Qc remains intact for tOr = 0. (The relation is valid 
even for 7^ upon the substitution w — > iiVr-) 
We examine the current conservation: 

i ^i°(T,x) = U{r)-H mr) - nN,f{^)] U{r) 



= [/(r)-i(-V,/(x) ^^^)u{r) 

= -Vfe?(T,x) . 



Hence 



^^i° + Vfc/ = • (19-39) 
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For the kernel V^" it implies that 



(19.40) 



Other relations such as T'^{p) = r^'*(— p) remain intact. The decomposition of 
r^'' is given by 



r"i'K,q) = gX 



r^^(t^,., q) = +iuJrqjIlQ - iejkQkTLf 



(19.41) 



^si^r, q) = -Sjk w^XIq' - ejkUrllf - {q''5jk - qjqk) 



.2rrB 



where all II^ 's are functions of and only. In a frame q = (g, 0) 



^ E 



(19.42) 



Mathematically 



r^^=r''^(ia;.,q) , U.'^ =Ilk{-ojt,q') . 



(19.43) 



As in the T = case, we need to evaluate one loop diagrams in Fig. 8 in Section 

10 to find . Computations arc completely parallel to those in Section 12. The 
diagram (a) yields the linear term in (19.32). The diagram (b) yields 



m 



m 



(19.44) 



For the diagrams (c), (d), and (e), the phase factor (j){x,y) in the propagator 
G{x,y), (19.27), completely cancells. We have 

rr = -^E/(^«o(p)Go(p-,) 

= E / {Go(p) ■ DjGoip -1) + DtGoip) ■ Goip - q)} 

=^^J(Sf {^^^'^p^ ^^^'^ - + ^^^^^ - 

+ D+D^Go{p) ■ Go{p -q) + Go{p) DtD+Go{p - q)} 

(19.45) 
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where 



/3 
2nr 

q = (Wr,q) , uJr = . 

Go{lOst'p) is defined in (19.29). Without confusion we have adopted the same nota- 
tion for the propagator as in the T = case. 

The only technical change to be made in comparison with the computations in 
Section 12 is the infinite sum over frequencies. Employing the formula 

CX) 

g{s) = — Res (tt cot nz ■ g{z) , aj) , 

s=—oo poles Oj of g(z) 

one easily finds that 



P Pn{l — pn) for Wr = and n = m, 

Pn Pm 



(19.46) 



otherwise. 



ILOr 



PniP) is the distribution function for the n-th Landau level. This is the only place 
where Pn shows up in the computation of T'^ . In other words, finite temperature 
effects in the linear approximation are contained solely in the discrete sum above. 
We shall see that the diagonal component at zero frequency, n = m and u>r- = 0, 
leads to unique behavior of anyon fluids at T 7^ 0. 

Working in the frame q = {q, 0), one finds H^'s to be 



g2n^ = r|° = -i-££/(a;.;n,m)CW(,z)2 , 



27r/2 

n — rn— 



gnf = -ir^2^^X^X^/(a;.;n,m)CW(gOCW(5Z) , (19.47) 

n=0 m=0 



^ 00 00 

n=0 m=0 

We are going to examine implications of the above result, particularly in the static 
case {u)r = 0), in the following three sections. 

20. Thermodynamic potential in inhomogeneous fields 

In this section we first summarize the result in the previous section in the form 
of free energy or thermodynamic potential for slowly varying static gauge fields 
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confugurations, from which it follows that the Meissner effect becomes partial at 
finite temperature at least in the self-consistent field method (SCF). Although it 
may be an artifact of the approximation which neglects vortices, the behavior found 
here is unique and seems essential for understanding properties of anyon fluids. 

The infinite sum in (19.47) can be performed for static, slowly varying gauge 
field configurations. Let us define 



oo 

=IZ(^+^rPn(l-Pn) • 
n=0 

It follows from (17.5) and (17.9) that 

Ao{T,B) = u 
Ao(T,0) = |7V| . 

At T = and B = 0, p„ = or 1 so that 

Ao = |iV| , Ai = — , A2 = ^^2 ' ••• 

Sn = , 

at T = , B = 

In evaluating off-diagonal sums (n ^ m) in (19.47), one also needs 

oo 

Y^In = S[I] 



n=0 



where 





In 






S[I] 




{Pn 


- Pn- 


_i)n 




-Ao 




{Pn 


- Pn- 


-i)n2 




- 2Ai 




{Pn 


- Pn- 


-i)n3 




-3A2- 


3A0 


[pn 


- Pn- 


_2)n(n- 


1) 


-4Ai 




[pn 


- Pn- 


-2) n(n - 


l)(2n- 1) 


- I2A2 - 


- 3A( 


[Pn 


- Pn- 


-3)n(n- 


l)(n-2) 


-9A2- 


fA 



(20.1) 



(20.2) 



(20.3) 



(20.4) 



With these preparations, we start to evaluate XIq' in (19.47) at Wr = 0. Employ- 
ing (19.46), one finds 

— R 
,^nf (0, e) = E ^^^('^O^ + 2^ E /^"(l - P^) C^M' (20.5) 

n^m n=0 
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where we have made use of e„ = (n + ^)/(ml^). For slowly varying configurations 
{ql <C 1) the expansions (12.30) - (12.32) can be employed: 



2-nErn 2\ \ ^ Pn Pm I r t i\2 



f^^ X n(n- 1) \ 4 

y-^On,m+l 0n,m+2J\ql) + ' 

OO f 



(gO + 5?; (90 



27r/2 ^''"v- 4 v^v ■ 32 

n=0 



so that 



nf and are similarly evaluated. One finds that for Ilf 



nf(o,a = +^{Ao-^Ai(50^ + . 



1 r 3 



For nf 

1 



g2nf(0,g2)=-- 



m 2-iTmP 



27rTO2/4 

so that 



(20.6) 



(20.7) 



Ao - 2Ai (qO' + {1^2 + I ^o) (<?0' + • • • } 



(20.8) 
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There are a few things to be recognized. IIq' develops a pole {l/q^) at T ^ 0, 
which, as we shall see shortly, leads to a partial Meissner effect. We also argue in 
the next section that it determines the scale of the phase transition temperature, Tc. 
Secondly, the diamagnetic term rie/m in g^nf' is cancelled by the induced term as 
at T = 0. No pole develops in Hf'. Thirdly Ilf (0, 0) determines the induced Chern- 
Simons term, which is not exactly N/2tt at T ^ 0. In other words, the cancellation 
between the bare and induced Chern-Simons terms is not exact. In some of the early 
literature in the anyon superconductivity it was said that the exact cancellation is 
essential for superconductivity, which, as we shall show, is rather misleading. 

It follows from (19.32) and (19.41) that 



fle[a, ^] = ^o[a, A]+ dx ne{ao + eAo) 




where llf{0,q'^) = Ilf{0,-W^). Insertion of (20.6) - (20.8) leads to 



f2e[a, A] = Oo[a,^] + dx < ne(ao + eAo) 






(20.9) 
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In applications at long wave length, dominant terms are given by 

=(const) + / I \{Fl + - ^ «o h^'^ 

-\iff CO («o + eAo)^ " " ci)(/o. + eF,,f (20.10) 



+ e(iV) ^ (1 - ci) (ao + eAo) (b^^) + eS^^)) 

ZTT 

where dimensionless constants Cj(T, B)'s are defined by 

47r2ne/3 27rne/3 T^rieP ,'on11^ 

Co = "^O , Cl = 5- dl , C2 = T— ()2 . (20.11) 

Equations are given by 

/¥y-o = :^ • (20-12) 

It is appropriate to examine numerical values of various coefficients. With the 
values m = 2me, rie = 2 x lO^^cm"^, and d = 5A, 

48 , — = 1.1 X 10-5 , ^ = 2.4 X 10-^ . (20.13) 



6p in c/j is approximately given by, for A'' = ±2, 



~ Max , e-^/-) ~ Max , 10-(--/-)) . (20.14) 

Other relevant coefficients are, for T = 100 K, 



= 1.2 X 10" , — |: = 14 . (20.15) 
2e^T ' 2mT ^ ^ 

All Cfe's at = are suppressed exponentially in the T — > limit, ci and 
C2 arc negligible (^ 1) for T < 200 K, whereas cq suddenly becomes large around 
T = 100 K. (See the discussion in Section 22, around (22.2).) The dominant finite 
temperature effect is contained in the cq term in (20.10), which represents an effect 
similar to the Debye screening in plasmas. 

21. Partial Meissner effect in SCF 

Anyon fluids have quite unusual behavior at T ^ 0. In this section we examine a 
response against static inhomogeneous external perturbations, both solving in real 
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configurations and looking at response functions. We write Eq. (20.12) in the 
form 

(Note that ak = —a'^.) The induced current -/^^(x) is given by 

jL(x)=e(^)|^(l-Ci)(6W+ei?(^)) 



— j CO (ao + eAo) - (1 - ci) ^/.(/ofc + eFofe) + ■ 



(21.2) 



_^(l_C2)e'='ai(6W+eijW)+--- . 
irm 

We remark that at finite temperature not only field strengths but also the time 
component of the vector potentials, ao + cAq, appears in the expression for in 
(21.2). 

The identities 

— e6(^) =div'E , 

% (21.3) 
^ e/o/c = dkB , 

may be employed to eliminate the Chern-Simons fields. Note that the integration 
of the latter leads to 

N 

— eao(x) = -B{x) + const , (21.4) 
27r 

where the constant has to be determined with the aid of, for instance, the neutrality 
condition at one point in a given configuration. (To be precise, only the constant 
part of ao + eAo is relevant in (21.2).) Substituting (21.3) and (21.4) into (21.2), 
one finds 

eJ,„, =(const) - (— j co^o + " (1 " ^O^^^J ^^^^ 



e4d={^(l-cr)-^(l-c.)}e-a,E 



(21.5) 



+ e(A^)-e'=' |(1 - c.)E, - (1 - c)-^ 9.(c^^.E) 

The equations in (21.5) correspond to the London equations in the conven- 
tional superconductors. Combined with the Maxwell equations, they determine 
electromagnetic fields inside anyon fluids. In general, however, one more equation, 
Eq. (18.2), has to be supplemented to fix the constant part of the magnetic field. 
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To illustrate the problem, we consider a configuration in which an anyon fluid 
occupies a half plane, say, xi > 0. We apply an uniform external magnetic field 
Bcxt in the empty space (xi < 0). The problem is to find -B(x) = B{xi) for a;i > 
with the boimdary condition B{0) — B^xt- 

To extract the essence, we suppose that ei?ext 6^"'. To good accuracy one 
can approximate v = \N\. With the aid of the numerical evaluation for various 
parameters given in Sections 18 and 20, one finds that the Maxwell equations become 



vi o Nth 
(1 + co)B - —„ — V^B - — — Co Aq + div E + (const) = 



(21.6) 



For the configuration under consideration one expects 

B(X) = Bin + (Bext - An) C-^^/^" {x^ > 0) . (21.7) 

Sin = B{+oo) does not vanish at T ^ 0. It is determined by Eq. (18.2) 

An = Sext+M(T,Bin) • (21.8) 

It follows from (21.6) that E2 = 0. 

It is checked posterior that \Ei/{B — -Bin)| <C 1 so that the first equation of 
(21.6) yields 

(1 + Co){B - Bin) - dl{B - Bin) = . 

Hence the damping length Aq is approximately given by 

Ao(T)^ = 4^ = -A1 (21.9) 

1 + Co e^rie 1 + Co 

and 

~-^~3x 10-^^(1 + co)V^ . (21.10) 
B{xi) - Bin toAo 

Due to the non- vanishing Bin, the damping length Ao(T) should not be confused 
with the penetration depth, which measures how fast the magnetic field decreases 
in the material. One complication in the calculation is that Eq. (21.8) cannot 
be solved analytically at T 7^ 0. We present the result of numerical evaluation in 
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Fig. 12 Magnetic field inside the anyon fluid in SCF 



As can be seen from the figure, Sin is vanishingly small at low temperature, but 
starts to increase around 70 K and becomes almost equal to Bext around 100 K. 
Here we have only a partial Meissner effect, at least in SCF. The magnetic field 
configuration is not a simple exponential decay. 

To avoid solving Eq. (21.8), one may apply a spacially alternating external 
magnetic field. In the linear response theory it is reduced to examinig Bext(a^i) = 
Bq €{xi) applied to a system occupying the whole space as we did in Section 16 at 
T = 0. We are going to show that the partial Meissner effect is observed in the 
response function, too. 

The equations to be solved are the same as in Section 16, with life's being 
replaced by n|''s. Employing the expansions (20.6), (20.7), and (20.8), and keeping 
dominant terms, one finds, instead of (16.9) and (16.10), 



As T changes, the response function for a charged anyon fluid smoothly varies. 
However, its behavior is different from that in conventional superconductors. 




(21.11) 



,22 



1 



1 + co + A g2 
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Recalling eJ^^^{ui,q) = —2Bq ■ (27r)^5(w)(5(<72), one finds 



■^tot ~ "^ext ^ "^ind ~ \ ^2 y -^ext 



t2 9 
CQ + A 2 
-2 - "^<=xt 



(21.12) 



1 + Co + A "g2 

The external current is not completely cancelled by the induced current: 

^t'ot(q = 0) = -r^Je'xt(q = 0)5^0 . (21.13) 

J- + Co 

The magnetic field is 



Bid) = f JL 
Qi 



_. f CO 1 , 1 1 1,2 

-'^ni + co??^l + co9? + Ao-M^^* 



(21.14) 



Co 9f + Aq ■ 

A new pole develops at q = 0. In the configuration space 

B(a;)=Boe(a;i)(--^ + -^6-1-^1/^4 . (21.15) 
[ i + Co 1 + Co J 

At T = 0, Co = so that the Meissner effect is complete. At T ^ 0, co 7^ 0, 
resulting a partial Meissner effect. As we have seen in the previous section, co(r) 
suddenly becomes very large around T = 100 K. Therefore the Meissner effect 
effectively terminates around this temperature. In the approximation (SCF) in use, 
however, there does not result a phase transition. We shall argue in the next section 
that a phase transition should result if vortices are incorporated. 

With the aid of (21.15) one can define an effective penetration depth, Ascf, 
which measures the rate of the change of the magnetic field. Ascf (I"; d) is related 
to the change of i3(a;i) over a distance d by 

e-dAscP = :^ . (21.16) 

■DO 

It depends on d. As a typical value we take d = \ . Then 

^ = -ln|-^ + -^e-(W/n . (21.17) 
Ascf ll + co 1 + co > 

Approximately 



Ascf rl + (e-l-5)co for co < 1 , 
A I Co for Co » 1 . 



(21.18) 
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The behavior of Ascf(T') is depicted in Fig. 13. 



Fig. 13 The penetration depth in SCF, AsCF defind in (21.17). The small tail in 
(A/Ascf)^ around 7=100 K is 
See the discussion in Section 22 



2 

(A/Ascf) around r=100 K is regarded as an artifact of the approximation. 



As we have demonstrated, the dominant finite temperature effect is contanied 
in co(T). The cancellation of the bare Chcrn-Simons term by the induced one, for 
instance, is not exact at T ^ 0, since ci(T) ^ 0. (See Eq. (20.10).) However, its 
effect is numerically negligible. For the Meissner effect, "cq" is important. It makes 
the Meissner effect partial at T 7^ 0. 

Before closing the section, we briefly mention about the subtlety in neutral anyon 
fluids. The response function for a neutral anyon fluid at flnite T is given by 



c-o + (m<zK) ^21.19) 
Co = e Co = do . 

It follows that the two limits, g — > and T — > do not commute with each other. 
It may reflect an instability in the system of neutral anyon fluids. 



22. Te 

How large is Tc, if there is a phase transition? Having analysed properties of 
charged anyon fluids, we are in an awkward position. In the linearized SCF, or 
equivalently in RPA, we have seen no evidence for a phase transition, or more pre- 
cisely, mathematical singuralities in physical quantities, at finite temperature. For 
instance, we have seen in the previous section that the penetration depth Xscf{T) 
rapidly increases around 100 K, but never diverges. 
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We argue that this is an artifact of the approximation in use, and that in a full 
theory a charged anyon fluid should exhibit a phase transition around T^. 

Crucially missing in the previous treatment is a vortex. It is missing, because the 
linearized version of the SCF equations (15.7) are linear in fields, and therefore do 
not admit a quantized flux. In terms of the effective theory obtained by integrating 
the fermion fields ip, one needs to retain higer order terms, cubic, quartic • • • in the 
field (a + e^)(i). It is a challenging problem to show how a vortex solution comes 
out from such an efi'ective theory. 

Previously the quantization of vorticity in a neutral anyon fluid was examined 
by Hanna, Laughlin, and Fetter in the Hartree-Fock approximation.^^ They showed 
that an elementary excitation has a vorticity given by (fundamental unit) /|A''|, 
although it has an infinitely large energy. Kitazawa and Murayama^^ have examined 
effects of vortex- antivortex pair formation in a neutral anyon fiuid at T 7^ 0. They 
have contended that there is a Kosterlitz-Thouless phase transition at T = | e for 
|A^| = 2. The underlying assumption is that there are vortex pair excitations with 
logarithmic interactions. 

Having vortex-antivortex pair excitatons is one promissing way of obtaining a 
phase transition in anyon fluids. Supposing abundant pair excitations, one still has 
to elaborate Kitazawa and Murayama's argument for charged anyon fluids. 

First of all interactions among vortices are not logarithmic at low temperature. 
The Meissner effect is operating so that interactions are exponentially suppressed 
at large distances. An energy of a single vortex, which is not known yet, must be a 
dominant factor at low T. 

The situation becomes more complicated as temperature increases. As we have 
observed in the previous section, the Meissner effect effectively terminates around 
T^. There would be no screaning of magnetic fields any more. The interaction among 
vortices becomes long-ranged, and the entropy factor becomes important. Whether 
or not this leads to a phase transition is a matter subject to hiturc investigation. 
Any way there is not any trace of superconductivity well above T^. It is quite likely 
that Tc, which separates the superconducting and normal states, turns out around 



signifies a temperature where co{T) becomes large. Prom (20.11) to (20.15) 

one finds 





(22.1) 



Numerically, for A'' 



±2, 



co(T) = 



121 K ni=^ _d_ 
T 2 X 10" cm-2 5A 




(22.2) 
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Typical values are 

T 50 75 100 125 150 175 200 

co(T) 2.4 X 10-6 1.6 X 10-2 1.2 15 81 260 610 

Note that if the values of both nf^^ and m are doubled, the value of cq is also 
doubled. 

The temperature dependence of cq is controled by the exponential factor. The 
critically important value is the ratio n'^r^ /m. With the given value in (22.2), cq 
suddenly becomes large around 120 K. 

With the assumption Tc ~ T'^ we conclude 

2 nf=2 2me ^ ^^qk . (22.3) 



\N\ 2 X 1014 cm-2 m 

Of course, the effective mass m is very difficult to determine experimentally so that 
one more piece of information is necessary to predict T^. One way is to express m 
in terms of the penetration depth at T = 0, A(0) = A , with the aid of (18.12) and 
(18.13). Then 

2 d (l^W,,^ . (22.4) 



" \N\^k \ A(0) 

We stress that the value Tc ^ 100 K is very natural in anyon fluids. The de- 
pendence Tc (X Ue or A(0)^2 has been observed in high Tc superconductors. This 
behavior, however, is not necessarily special to anyon superconductors. 

Although the discussion in this section is only plausible and further investigation 
is necessary, one might take (22.4) as a very encouraging result. 

23. Other important issues 

In this article we have analysed some of the basic problems in anyon fluids, at- 
tempting to summarize the first two years of the theory of anyon superconductivity. 
We have seen that various approaches are equivalent, leading to many interesting 
physics consequences. Nevertheless, it is appropriate to say that so far we have had 
only partial understandings of the full theory. 

There are many important issues left over. We list them for readers' convenience. 
For details readers should consult original papers and other review articles. 

23.1. Beyond RPA and the linearized SCF 

Going beyond RPA and the linearized SCF is important in many respects. RPA 
and the linearized SCF fail to predict a phase transition. It is essential to incorporate 
vortices in the theory. 

Hanna, Laughlin, and Fetter^** examined various quantities in the Hartree-Fock 
approximation in the first quantized theory, as was described in Section 7. They 
have found that, for | A''] =2, the correction to the phonon spectrum in neutral anyon 
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fluids in the long wave length limit is relatively small (~ 10 %). More recently 
Dai, Levy, Fetter, Hanna, and Laughlin^^ have performed a diagram analysis at 
r = equivalent to the Hartree-Fock approximation. In addition to recovering 
the old result, they have shown that there results an important modification in 
the behavior of the spectrum at short wave lengths. Furthermore Cxy vanishes for 
TV = ±1 thanks to higher order radiative corrections. RPA and the linearized SCF 
predict a non- vanishing axy even for = ±1, which is certainly wrong as fermions 
are merely converted to bosons, but not to genuine anyons. 

Dai et al. have solved the Schwinger-Dyson equations numerically which in- 
volve about 100 different diagrams. They have noticed the importance of gauge 
invariance, and have observed many cancellations among various diagrams. Their 
Feynman rules are based on the Hamiltonian obtained after eliminating Chern- 
Simons fields. As we have recognized in Sections 9 - 11, keeping the Chern-Simons 
gauge fields as auxiliary fields greatly simplifies computations. Gauge invariance 
is easily implemented, and the notion of self-consistent fields (beyond the linear 
approximation) can be established. 

23.2. Pair- correlation 

Is there "Cooper" pairing in anyon superconductors? The Hartree-Fock ground 
state employed in the literature (see Sections 6 and 7) does not look like the BCS 
ground state. Is there a different kind of pairing, then? Is there an off-diagonal 
long range order? The answer has not been known for sure. 

There is an indication for pairing in the A'' = ±2 theory, which, however, is quite 
different from the Cooper pairing. The unique feature of the Hartree-Fock ground 
state is that the complete filling of the Landau levels is achieved independent of the 
density rig, provided that N is an integer (^ 0). 

To be precise, suppose that each Landau level has Nl available states so that 
the total particle number is = Ug ■ (vol) = \N\ ■ Nl- Assume that |A^| > 2. If 
one tries to add or delete one particle to or from the system, one necessarily has 
to put the particle in the next level, or make a hole in the top filled level, in order 
to preserve the Landau level picture. In other words, the picture of the complete 
filling breaks down. 

However, if a set of |A^| particles are added or deleted, one can still maintain the 
complete filling. One of the Chern-Simons field equations, {N/2Tr) b = j'^, implies 
that the increase (decrease) of the particle number leads to the increase (decrease) 
of the Chern-Simons magnetic field such that precicely one more (less) state is 
available in each Landau level. 

The states with the particle number N^, and ± lA^] are very much alike. In 
a macroscopic system A^g ^ 1, thermal fluctuations give AA'e ^ VNe- It is quite 
likely that the real ground state is not an eigenstate of the particle number, but is 
a coherent state: 

*g(^)= Y1 e"'^ ^G{Ne + k\N\) . (23.1) 

k 
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In particular, for N = ±2, the structure (23.1) is exactly the same as in the BCS 
theroy.^°^ 

It is not clear if (23.1) implies pairing in the N = ±2 theory. We should 
remember, however, that the structure of the coherent state is indispensable in 
understanding many phenomena in superconductivity. 

23.3. Flux quantization 

A magnetic flux is trapped by a superconducting ring. The flux takes quan- 
tized values in the unit of 27rS.c/2e. Can the N = ±2 anyon theory explain this 
behavior? No convincing argument has been provided. Leggett^^ has argued that 
the flux quantization is not achieved in anyon theory at least in the Hartree-Fock 
approximation. 

One needs to show two things. First it must be shown that an energy is locally 
minimized when a flux takes a quantized value. Secondly, the energy barrier height 
between the states with no flux and with one unit of flux is proportional to the 
volume, but not the boundary area, of the superconducting ring. 

23.4. Vortices 

We have often mentioned in the preceeding sections that the establishment of 
vortices in anyon superconductors is one of the major problems to be solved. The 
previous analysis by Fetter, Hanna, and Laughlin must be elaborated. Inclusion 
of electromagnetic interactions is essential to take account of the Meissner effect. 
Nonlinear terms in SCF must play an important role. At finite temperture vortices 
might lead to a phase transition, too. 

23.5. The Josephson effect 

The nature of the coherent state (23.1) is most important for the Josephson 
effect. A Josephson junction consists of two superconductors separated by a bar- 
rier. Electron tunneling through the barrier brings about phase coherence over the 
entire system. The energy is minimized by volume if the two phases, which char- 
acterize 9's in (23.1) of the two superconductors, coincide. The difference between 
the two phases should generate a current. 

Experiments show that in order for anyon superconductivity to describe high Tc 
superconductors, TV must be equal to ±2. It seems that a Josephson effect should 
exist even for a junction between a BCS superconductor and an anyon supercon- 
ductor of TV = ±2. 

23.6. Interlayer couplings 

High Tc superconductors have the layered structure characterized by CuO 
planes. Our analysis has been performed in the effective two-dimensional theory ob- 
tained by the dimensional reduction. The implicit assumption was that the system 
is uniform in the direction perpendicular to the two-dimensional CuO planes. 

Anyon theory has a parameter N, or the generated statistics phase ^statistics = 
-k/N. Physics depends on exp(i^statistics)- The theory with A'' differs from that with 
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—N if I A'' I > 2. The two theories are related by P (parity) and T (time reversal) 
transformations. The idea behind the degeneracy of the ground state is that P and 

T symmetry is spontaneously broken. 

P and T invariant quantities such as the penetration depth and resistance do 
not depend on the sign of N. Moreover we have seen that even though the Chern- 

Simons magnetic field b/e is very large 1000 T), the dependence of the P- and 
T-odd magnetization M{B) on the Maxwell magnetic field B is symmetric to good 
accuracy, M{—B) ~ M{B). The asymmetry arises only to the order 0(eB/b). (See 
Section 18.) 

Nevertheless it is important to know how the sign of is ordered among adja- 
cent layers. Is it ordered ferromagnetically with the same sign (FM ordering), or 
antiferromagnetically with the alternate sign (AFM ordering)? Or, is it randomly 
distributed? 

This problem has been examined by Rojo, Canright, and Leggett.^°^ Interactions 
among electrons in different layers fix a pattern of the ordering. There are two types 
of interactions. One is of a potential type, and the other is the hopping of electrons 
from one layer to adjacent ones. 

The detailed examination in the case of potential interactions has been provided 
by the above authors both numerically and analytically. They have shown that 
T-invariant potential interactions always prefer the AFM ordering. The hopping 
interaction is expected to induce a Josephson efi'ect and lead to the FM ordering. 
It is not clear which one is dominant. 

23.7. P and T violation 

Many experiments have been performed to check P or T violation in high Tg 
superconductors. The result is confusing, but a fair statement is that so far there 
has been no solid evidence for P or T violation in high Tc siipcrconductors. 

As explained above, P and T are ordered either ferromagnetically or antiferro- 
magnetically among layers. Most of the experiments done so far measure P and T 
violation in bulk. Therefore, if material has the AFM (P-, T-) ordering, the effect 
cancells in bulk. 

Experiments in this category are the electromagnetic wave polarization and Hall 
voltage. Originally proposed by Wen and Zee, the polarization experiment tries to 
measure the P, T violation effect, detcrminig the transmission and reflection co- 
efficients of injected polarized electromagnetic waves. The polarization vec- 
tor is rotated in anyon superconductors. There is inconsistency among various 
experiments, however. We note that even though Wen and Zee argue that there 
must be appreciable effects in FM-ordered anyon superconductors, a really mi- 
croscopic computation of the magnitude of the effect is still lacking. Also some 
criticisms have been provided on the interpretation of experimental results. 

The Hall voltage experiment measures the temperature dependence of the trans- 
verse voltage (Hall voltage) when a current flows in a thin fllm.^^'^°^'^^° There is 
a microscopic calculation of the effect. Theory predicts a peak in the Hall voltage 
around Tc- The peak value predicted is, for a thin film of thickness 1000 Awith a 
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current 10~^ Amp, Vua,n ~ 2 x 10"'' Volt. It is inversely proportional to the thick- 
ness. The effect is tiny. Preliminary experiments have been performed. Due to the 

inhomogeneity of samples and also tiny temperature variation in the samples, no 
conclusion has been obtained concerning the existence or non-existence of the Hall 
voltage."^'^^^ 

There is one experiment which measures P, T violation in one layer, and there- 
fore is sensitive even for AFM-ordered anyon superconductors. It is the muon spin 
relaxation experiment. Injected muons are stopped in high Tc superconductors. 
Since muons are charged, the distribution of electrons is deformed. In effect, the dis- 
tribution of holons, or our xj) particles, deviates from the uniform value. (5J°(x) ^ 
results. In anyon superconductors it induces a current J^^j ^ 0, since ^ 0. 
Jj^ J ^ in turn generates Maxwell magnetic field Bf^^, which is felt by muon spins. 
Muon spins start to preccss, which can be observed experimentally. 

Halperin, March-Russel, and Wilczek gave a plausible argument, predicting 
And ~ ^- '^'^^ experiment observed no effect. We remark that the mag- 
nitude of the effect can be determined more microscopically from the knowledge of 
the response function Qc at both T = and T without making the long wave 
length approximation. 

23.8. Anyons in spin systems 

In this article we have not discussed how anyon excitations arise in material, 
particularly in spin systems. We have started with the picture that there are exci- 
tations called "holons" which obey half-fermion {N = ±2) statistics. 

The derivation of anyons, or fractional statistics, in realistic spin models for 
high Tc material has been attempted by many authors. The issue has not been 
settled yet, since the arguments involve many approximations for which justification 
is not clear. Readers are advised to read original papers. A closely related subject 
is the existence of the flux phase or chiral spin liquid state. 

We note that Laughlin has given a spin model which has Laughlin's wave func- 
tion in the fractional quantum Hall effect as an exact ground state wave function. 
This model provides the exsistence proof of anyons in spin models. 

23.9. Variations of anyon models 

We have analysed one particular anyon model, namely non-relativistic spinless 
fermions with the minimal Chern-Simons interaction. It is the simplest model of 
anyons, and is based on the holon picture of Anderson's. 

There are many variations. They are interesting in their owm right. Historically 
Chern-Simons gauge theory was first analysed in relativistic field theory. Along 
this tradition Lykken, Sonnenschein, and Weiss have examined a relativistic anyon 
model, Dirac fields with Chern-Simons interactions. ^^'^^ They have argued that 
the neutral model retains a superfluidity to all orders at T = 0. In passing, Imai 
et al. have shown that the non-renormalization theorem for the induced Chern- 
Simons coefficient holds even in non-relativistic theory. At finite temperature the 
relativistic model behaves in a fashion similar to, but not same as, the nonrelativistic 
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model. To discuss a Meissner effect, superconductivity etc. in condensed matter 
systems, one has to analyse nonrelativistic models. Based on non-vanishing finite 
temperature corrections to the induced Chern-Simons coefficient. Lykkcn et al. have 
incorrectly concluded that a superconductivity is lost at T 7^ 0. As we have seen in 
Sections 20 and 21, the important finite temperature correction is the co term, but 
not the ci term (Chern-Simons cocfHcicnt). in the non-relativistic theory. 

We have supposed that particles (anyons) have a single component, i.e. they 
have the same coupling to the Chern-Simons fields. There might be two kinds of 
anyons, a half of them having the + coupling and the other half having the — 
coupling. Furthermore, in addition to the minimal gauge coupling to Chern-Simons 
fields, particles might have magnetic moment interactions. Such a model has been 
investigated.^'^ 

So far we have started with fermion fields ip. ft is also possible to start with 
boson fields. It is not exactly the same as the fermion model, since a bose field 
can condensate by itself and there is no complete filling of Landau levels. As was 
pointed out by Boyanovsky et al., boson models have rich structiires many of which 
need to be clarified further. It is also known that boson models are particularly 
useful to construct phenomenological theory of fractional quantum Hall effect. 
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